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SUMMARY 


The  principal  objective  of  the  Flameholder  Combustion  Instability  Study  was  to  develop  an 
analytical  model  of  the  flame  stabilization  and  propagation  process  in  a conventional  turbofan 
engine  augmentor  and  describe  the  response  to  low  frequency  pressure  )scillations.  The  model  is 
to  be  used  for  identification  of  the  driving  mechanisms  of  low  frequency  instability  (rumble)  in 
turbofan  augmentors. 

The  program  was  conducted  in  two  phi._,es.  Phase  I represented  the  bulk  of  the  activity  and 
encompassed  development  of  the  analytical  model  and  execution  of  subscale  experimental  rig 
testing.  Phase  11  included  compute. ization  of  the  model  in  a format  suitable  for 
analysis/prediction  of  a turbofan  augmentor  and  comparison  of  the  results  with  the  available 
liteiature  and  engine  data.  The  concepts  were  also  extended  to  include  two  novel  augmentor 
types  with  prediction  of  their  relative  stability. 

The  Phase  I model  relied  on  analysis  of  the  individual  processes  which  control  the  flame 
stability  in  low  inlet  temperature  airstreams.  Since  the  fuel  is  highly  liquid,  a model  was 
conceived  for  stabilization  in  a two-phase  fuel  spray  where  the  flameholder  serves  the  dual  role 
of  controlling  both  vapor  composition  and  kinetic  stability  limits.  The  response  is  quasi-steady 
and  represents  the  change  in  augmentor  efficiency  with  fuel-air  ratio  and  inlet  conditions.  The 
model  evaluates  the  various  interrelated  effects  of  injector,  flameholder  and  augmentor  design. 
The  results  of  the  experimental  study  were  integrated  mto  the  analyses.  The  model  predictions 
are  unique  in  terms  of  rich  limit  predictions  and  flow  variable  response  factors. 

The  concept  was  computerized  to  analyze  a turbofan  engine  and  agreed  very  well  with  test 
experience  on  tliat  engine.  The  behavioral  characteristics  were  in  concert  with  available  data  and 
other  models. 

The  two  novel  augmentors  were  found  to  provide  increased  stability  through  a reduction  in 
reliance  on  fuel-air  ratio  for  flame  speed,  removal  of  the  liquid  vaporization  dependence,  and 
increased  dynamic  stiffness  of  the  piloting  source  relative  to  flameholders. 
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SECTION  I 


INTRODUCTION 


'Php  phrnttnu'non  of  romhusfion  inxlahility  has  plnpirH  all  ly|K's  of  hinh  heat  n-lcaM' 
ct'mhuslors  from  industrial  furnnct's  to  nn'ki't  rnjtinps.  In  ftpnf'ral,  howpvpr,  th<’  prohlnn  has  hern 
most  srvorp  in  fli^'ht  propulsion  systmu.  such  as  turlaijrts,  ramjets,  .and  r<K'kcts  where  weight 
ronsideraf ions  dictate  hichly  efficient  structures.  In  such  applications,  the  preK.surc',  vihralion, 
and  heat  loads  resulting:  from  lomhustion  instability.  sufK’rimjKised  on  the  normal  loading,  are 
usually  dc.-.tructivp. 

In  airhreathing  engine.!,  high  frerpiency  rombustion  instability  problems  were  First 
encountered  in  the  early  I9.'i()'s,  and  solutions  were  sought  through  mathematiial  modeling  and 
analytical  studies  directed  toward  an  understanding  of  the  phenomenon.  I'nfortunateiv.  the 
computer  technology  and  enalytical  techniques  of  20  years  ago  |)roved  inadequate,  and  “cut  and 
try"  empirical  approaches  involving  changes  in  flameholders,  ci'mbiistion  chamber  shafie,  fuel 
injectio.-,  velocity  jirofiles,  and  flame  piloting  wen’  attempted.  Fuel  additives  and  combustion 
rhcmtier  t affies  were  also  tested.  Allhoug!i  some  of  these  apfiroaches,  notably  baffles,  jirorluced 
marginal  improvement,  the  problem  was  not  solved  until  danqiing  devices  in  the  form  of 
acoustical  absorbers  (screech  liners)  were  introduced.  Screech  liners  are  used  routinely,  and  high 
frequency  instability  is  no  longer  regarded  as  a problem. 

Larger  and  more  powerful  turbopropulsion  systems  are  presently  Ireing  designed  and 
developed.  Hecause  of  the  large  physical  dimensions  of  augmentors  used  in  these  systems,  their 
natural  acoustic  motles  have  corresjamdingly  long  wavelengths;  therefore,  combustion  instability 
ran  wcur  at  very  low  frequencies,  i.e.,  approximately  200  Hz  or  less. 

The  oc  urrence  i.f  instability  at  lower  frequencies  makes  use  of  screech  liners  of 
conventional  design  difficult.  To  obtain  adequate  damping,  the  absorbing  devices  are  designed  so 
that  the  resonant  frequency  corresponds  to  tlie  frequency  of  the  expected  mode  of  instability.  The 
required  cavity  volume  is  inversely  pro|)ortiona!  to  the  square  of  the  resonant  frequency; 
therefore,  low  frequencies  require  large  volumes.  Large  cavity  volumes  can  be  accommodf.ted  by 
increasing  the  augmentor  envelope,  but  this  produces  an  unacceptable  increase  in  engine  weight. 
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With  the  introduction  of  military  augmented  turbofan  engines.  P&WA  became  involved  in 
the  prihlems  of  very  low  frequency  instability.  The  problem  of  these  very  low  frequencies,  called 
rumble,  have  lieen  reduced  through  combined  experimental  and  analytical  techniques.  This 
experience  has  emphasized  the  necessity  to  understand  the  fundamental  mechanisms  involved  in 
order  to  formulate  a meaningful  analytical  effort  and  the  necessity  to  relate  this  effort  to  physical 
hardware  and  processes. 

Mathematical  mixlels  and  analytical  studies  of  low  frequency  instability  have  received  less 
attention  than  studies  of  high  frequency  instability  and,  unfortunately,  have  not  been  notably 
successful.  The  problems  of  low  frequency  instability  have  been  solved  or  circumvented  by 
empirical  methods.  For  rumble,  the  development  problem  is  even  more  expensive  and  time 
consuming  than  for  screech  since  rumble  usually  occurs  only  at  high  altitude,  low  flight  Mach 
number  operation.  This  requires  the  extensive  use  of  altitude  simulation  teat  facilities  to  develop 
a stable  augmentor  for  a turbofan  engine. 

The  main  feature  which  limits  the  usefulness  of  many  combustion  stability  models  is  the 
absence  of  a direct  correlation  between  the  physical  hardware  of  n real  augmentor  and 
combustion  stability.  Experience  has  shown  that  relatively  subtle  r.lterations  in  flameholder 
geometry  ern  produce  profound  changes  in  rumble  limits.  For  a usable  design  and  evaluation 
tool,  the  model  must  he  able  to  relate  directly  to  such  geometry  changes. 
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Analyli'  '1  and  experimental  studies  we'e  conducted  to  formulaic  the  model  necessary  to 
determine  the  response  of  an  auRmentor  flameholder  to  velocity  and  pressure  oscillations  typical 
of  low  frequency  instabilities.  The  influences  evaluated  include; 

1.  Approach  velocity,  pressure,  temperature,  and  turbulence 

2.  Flameholder  geometry  and  blockage 

3.  Fuel  distribution  and  vaporization 

4.  Heat  addition  to  flameholder  wake 

6.  Inclined  flameholder. 

Special  attention  was  given  to  the  processes  of  dynamic  flame  stabilization  in  cold 
airstreams  with  partially  vaporized  liquid  fuels.  This  area  is  one  which  is  unique  to  turbotan 
engines  and  one  which  has  demonstrated  the  largest  influence  on  the  occurrence  and  severity  of 
rumble. 

The  experimental  studies  of  the  flame  stabilization  process  and  the  dynamic  pressure 
response  of  the  flameholder  and  augmentor  system  were  conducted  using  the  existing  PC:WA 
Rumble  Simulator  Rig.  This  test  rig  was  specifically  designed  to  investigate  the  mechanisms 
involved  in  rumble. 

1.  PHASE  I ANALYTICAL  APPROACH 

The  approach  utilized  during  Phase  I of  this  program  was  to  formulate  a model  for  the 
process  of  flame  .stabilization  and  combustion  in  a turbofan  augmentor  in  terms  of  basic  physical 
phenomena  related  to  the  aerodynamics  and  thermodynamics  of  the  fuel  preparation  and 
combustion  processes.  By  basing  the  model  on  evaluations  of  individual  processes  which  cause 
the  rumble  effect,  rather  than  correlating  the  effect  to  the  overall  engine  parameters,  several 
significant  improvements  are  available: 

• Increased  insight  on  the  fundamental  causes  becomes  available. 

• The  individucl  analysis  may  draw  from  an  experimental  data  base  which 
fully  spans  the  ranee  of  operating  parameters  where  rumble  occurs. 

The  model  which  has  been  formulated  is  a steady-state  analysis  of  the  processes  which 
control  the  stabilization  and  combustion  in  the  turbofan  augmentor.  Complete  details  are 
presented  in  Sect  ion  II  of  this  report.  A summary  is  presented  below  as  an  outline  of  the  solution 
procedure.  The  model  introduces  a concept  of  stabilization  in  the  fan  duct  section  of  the 
augmentor  where  the  process  of  liquid  film  vaporization  from  the  flameholdtr  surface  is  a critical 
item  in  controlling  the  recirculation  zone  vapor  phase  fuel-air  ratio. 

The  model  analyzes  a two-dimensional  streamtube  (Figure  1)  and  all  comments  are  related 
to  this  streamtube.  Since  different  approaches  are  taken  for  the  fan  and  core  streamtubes,  they 
are  discu8,sed  separately. 

a.  Fan  Duct  Streamtubes 

The  effect  of  inlet  air  temperatures  which  are  significantly  lower  than  turbine  exhaust 
temperatures  is  to  reduce  the  liquid  fuel  vaporization  level  to  a degree  where  these  processes  must 
be  closely  studied.  The  degree  of  vaporization  of  the  droplets  of  fuel  is  so  low  between  the  injector 
and  the  flameholder  that  the  vapor-phase  fuel-air  ratio  is  below  the  lean  flammability  limit. 
Since  only  gaseous  species  may  traverse  the  shear  layers  and  enter  the  recirculation  zone  wake  of 
the  flameholder,  there  would  not  be  sufficient  fuel  available  to  provide  a stable  flame  if  only 
droplet  vaporization  were  present. 
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’riu'  liiodol  proposes  that  liquid  fuel  droplets  impart  the  surface  of  the  flameholder  in 
suifiririit  (piantity  to  form  a surface  fdin  of  liquid  phase  fuel.  This  thesis  is  substantiated  by 
measurements  of  tiie  ilameholder  surface  temperature  in  a turbofan  engine  which  revealed 
signilicant  decreases  when  ti  c fuel  was  rejected.  The  thesis  further  supposes  that  this  liquid  film 
undergoes  |iartial  vnporizat io.a  due  to  heat  transfer  from  the  reacting  wake. 

'riii.s  partin’  svirface  vaporization  provides  the  necessary'  source  of  vapor-phase  fuel  to 
provide  a llammahle  mixture  in  the  wake. 

The  wake  recirculation  rate  is  evaluated  from  empirical  data.  The  reaction  rate  within  the 
wak('  is  coaitled  with  the  recirculation  rale  through  a stirred  reactor  analogy’  to  evaluate  the  wake 
reaction  efficiency  and  temperature. 

Downstream  of  the  wake  region,  a turbulent  flame  is  initiated  in  the  shear  layers  and  its 
propagation  rate  calculated.  The  tlame  speed  is  related  to  the  degr  ■ of  surface  vaporization  and 
wake  react  ion  to  account  for  the  local  quenching  caused  by  nonvaporir  d fuel  traversing  the  shear 
layer  zone. 

The  ultimate  degree  of  tlame  penetration  into  the  approach  mixture  determines  the 
streamtube  efficiency.  Multiple  streamtubes  are  mass  averaged  to  evaluate  overall  efficiency. 

Tile  model  consists  of  .separate  analyses  of  each  of  these  processes.  Details  are  presented  in 
Section  II. 

b.  Coro  Streamtubes 

The  much  higher  turbine  exhaust  gases  provide  complete  droplet  vaporization.  Thus,  the 
wake  fuel-ttir  ratio  is  assumed  to  be  the  same  as  the  total  injected  fuel-air  ratio.  As  before,  the 
wake  recirculation  rate  is  evaluated  from  empirical  data.  The  same  stirred  reactor  analysis  is 
used  to  evaluate  the  wake  reaction  efficiency.  Turbulent  flame  penetration  into  the  approach 
mixture  is  calculated  to  evaluate  the  streamtube  efficiency. 

2,  Eiv.-Zni.MENTAL  APPROACH 

The  experimental  approach  selected  was  a systematic  evaluation  of  the  influence  of  various 
geometric  and  flow  variables  under  controlled  conditions  which  simulated  conventional 
augmentor  operation  at  high  altitude  and  low  flight  Mach  number.  The  test  apparatus  used  for 
this  study  was  the  Pratt  & Whitney  Aircraft  Rumble  Simulator  Rig.  This  test  apparatus  is 
cs.seiitially  a long  duct  of  12  in.  ID  which  contains  provisions  for  installing  flameholder  units, 
spraybars  and  viewing  ports.  It  is  modularized  to  allow  alterations  in  upstream  length  and 
combustion  length.  It  includes  an  upstream  choked  orifice  plate  to  simulate  the  fan  or  turbine 
exit  plane  and  a variable  area  exhaust  nozzle. 

Tlic  rig  is  fed  airflow  from  the  bleed  on  a slave  JT-4  engine  and  exhausts  through  an  ejector 
.system.  Nonvitiated  inlet  air  is  available  to  temperatures  over  600°F  and  vitiated  air  to 
temperatures  over  1500°F.  Subatmospheric  test  chamber  pressures  are  possible.  A full 
descri[)tion  of  the  test  apparatus  is  given  in  Section  II,  2.  of  this  report. 

'I  he  experimental  program  was  run  to  evaluate  the  relative  influence  of  the  following 
parameters  on  the  combustion  efficiency  and  rumble  intensity: 

3 Fuel-air  ratio 
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• Flameholder  geometry 

— Blockage  and  width 
— Gutter  apex  angle 

• Flameholder  to  spraybar  separation  distance 

• Drafted  flameholder  gutters 

• Upstream  (duct)  length 

• Combustion  length 

• Flow  variables 

— Static  pressure 

— Approach  velocity  and  temperature 
— Free-stream  turbulence  intensity 
— Fuel  vaporization 
— Heat  addition  to  flameholder  wake 
— Fuel  distribution 

These  tests  were  run  in  two  series  over  a six  week  period.  Fourteen  separate  configurations 
were  run  over  a flowrate,  pressure,  temperature  and  fuel-air  ratio  excursion  test  matrix.  Extensive 
static  and  dynamic  instrumentation  was  utilized  as  well  as  a senes  of  high  speed  motion  pictures. 

The  basic  objective  of  this  approach  was  to  reduce  the  amount  of  uncontrolled  variables  to 
a minimum  and  try  to  limit  test  changes  to  one  variable  at  a time. 

3.  PHASE  I RESULTS 

The  model  for  flame  stabilization  and  propagation  in  turbofan  augmentors  was  formulated. 
The  analytical  predictions  matched  the  results  of  the  experimental  program  well.  The  model 
correctly  analyzed  the  impact  of  two-phase  fuel  on  the  combustion  behavior  of  the  system.  The 
combined  analytical  and  experimental  programs  produced  the  following  results: 

• Augmentor  rumble  may  be  treated  as  a classical  combustion  stability 
problem,  where  the  oscillation  is  driven  by  the  pressure  or  velocity  response 
of  the  overall  augmentor  heat  release  rate. 

• The  overall  heat  release  rate  is  sensitive  to  variations  in  velocity  and  pressure 
not  only  through  the  flame  spreading  in  the  free-stream  but  also  through  the 
influence  of  the  flameholder  wake  region  on  the  downstream  flame  speed. 

• The  ultimate  level  of  augmentor  efficency  is  very  sensitive  to  the  nature  of  the 
composition  and  combustion  processes  in  the  flameholder  near  wake  region. 

• The  flameholder  wake  region  composition,  i.e.,  fuel-air  ratio,  is  not 
necessarily  the  same  as  the  overall  fuel -air  ratio  and  is  very  sensitive  to  the 
exact  geometry  of  the  flameholder  as  well  as  operating  aerothermodynamic 
conditions. 
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• The  combustion  model  correctly  predicted  the  influences  of  the  major 
variables  on  the  overall  combustion  efficiency,  e.g.: 

— Total  fuel-air  ratio 
— Flameholder  approach  conditions 
— Vaporization  level 
— Flameholder  blockage  and  width 
— Wake  heat  addition 
— Turbulence  level 

• The  cause  of  the  increased  rumble  sensitivity  with  increases  in  fuel-air  ratio 
and  approach  flow  severity  was  identified  as  a decrease  in  wake  reaction 
efficiency  and  flame  speed. 

• The  following  will,  in  general,  reduce  rumble  (with  the  converse  true); 

— Wake  heat  addition 

— Decreased  wake  fuel-air  ratio  through  fiameholder  geometry 
— Increased  turbulence  level 
— Increased  fuel  vaporization 

• In  addition,  the  following  produce  more  severe  rumble: 

— Increased  spraybar  to  flameholder  separation 
— Increased  fuel-air  ratio  or  aerodynamic  loading 

4.  PHASE  II  APPROACH  AND  RESULTS 

The  two  tasks  performed  during  Phase  11  of  this  contract  consisted  of  computerization  and 
evaluation  of  the  analytical  model  for  combustion  in  a conventional  mixed  flow  turbofan 
augmentor  and  preliminary  extension  of  these  concepts  to  alternative  augmentation  concepts. 
Along  with  this  latter  task,  predictions  of  the  stability  of  these  alternative  concepts  was 
performed. 

The  computerization  of  the  analysis  from  the  Phase  I formulations  proceeded  in  a relatively 
straightforvv-'''d,  if  somewhat  time  consuming,  manner.  The  computer  program  utilizes  a 
supervisory  main  with  the  analysis  performed  in  a subroutine  set.  This  approach  allows  rapid 
alteration  and  improvement  on  the  individual  process  calculations  without  disturbing  the  main 
program  logic. 

The  computer  program  will  analyze  the  flame  stabilization  and  propagation  process  behind 
a bluff  body  stabilizer  in  either  the  fan  duct  or  core  stream  zones  of  a mixed  flow  turbofan 
augmentor.  Different  sequences  ot  subroutines  are  exercised  for  these  two  cases.  The  lower 
temperature  fan  duct  analysis  performs  a two-phase  fuel  stabilization  calculation  while  the  high 
temperature  core  stream  case  is  analyzed  on  the  basis  of  fully  vaporized  fuel. 

The  program  has  been  exercised  to  evaluate  its  predictions  against  the  results  of  the 
experimental  program  and  historical  data.  Predictions  have  been  made  for  a turbofan  engine 
augmentor  and  compared  with  test  experience  on  that  engine  development  program.  The  results 
of  these  studies  are  very  encouraging  in  the  following  manner: 

• The  behavior  of  the  efficiency  versus  total  fuel-air  ratio  exhibits  the  observed 
sharp  rise  from  the  lean  limit,  peak  on  the  lean  side  and  slower  decline  until 
rich  blowout. 
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• The  influence  of  wake  heat  addition  is  correct  with  respect  to  efficiency 
increase  magnitude  and  rich  limit  increase  when  compared  to  either  test  data 
OT  engine  experience. 

• The  model  pr>.  diets  augmentor  blowout  very  well  when  run  for  the  test  case 
engine  and  compared  with  development  experience. 

• The  quantitative  values  and  the  variation  with  altitude  are  in  agreement  with 
experience. 

• The  influence  of  inlet  variables  is  correct  in  direction  and  approximate 
magnitude. 

The  efforts  under  Task  II  of  Phase  II  were  performed  on  extension  to  alternate  concepts  of 
turbofan  augmentor  design  and  evaluation  of  the  possible  rumble  stability  improvement  due  to 
these  concepts. 

Rumble  has  been  identified  as  the  result  of  the  dynamic  response  of  the  augmentor 
e^ciency  to  variations  in  (in  order  of  sensitivity): 

• fuel-air  ratio 

• pressure 

• velocity 

• temperature. 

For  such  a situation,  rumble  may  be  reduced  by  any  design  action  which  reduces  the 
response  of  the  dynamic  efficiency  variation  to  any  of  these  input  variables.  Such  efficiency  is 
essentially  a measure  of  the  rate  cf  transverse  flame  propagation,  rumble  is  reduced  by  any  design 
or  concept  variation  which  reduces  the  dependence  of  flame  speed  on  these  four  variables. 

Two  alternate  concepts  of  augmentor  design  have  been  advanced  in  recent  years.  These  are 
the  Swirl  Augmentor  and  the  VORBIX  Augmentor.  Both  of  these  attempt  to  provide  forced 
mixing  mechanisms  which  reduce  the  response  of  efficiency  to  fuel-air  ratio  dynamic  variations. 
The  effect  of  forced  mixing  is  to  remove  the  sensitivity  of  transverse  flame  speed  to  the  local  heat 
release  rate  (i.e.,  fuel-air  ratio)  and  replace  it  with  a turbulent  mixing  function  which  is  primarily 
a function  of  the  design. 

Such  in  approach  will  produce  an  augmentor  design  which  is  less  prone  to  rumble.  The 
problems  which  still  exist  in  these  concepts,  however,  are  the  piloting  response  and  fuel 
distribution.  In  general,  severe  pilot  response  to  flow  variations  or  severe  nonuniformities  in  fuel 
distribution  will  produce  higher  rumble  response  characteristics.  These  topics  are  treated  in 
Section  II. 
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TECHNICAL  DISCUSSION 
1.  PHASE  I MODEL  FORMULATION 
a.  Rumble  Mechanism 

During  Phase  I a combustion  analysis  was  generate  for  the  process  of  stabilization  and 
propagation  of  flame  in  the  presence  of  liquid  hydrocarbon  fuels  with  conventional  bluff  body 
stabilizers.  The  onset  of  rumble  has  been  identified  as  occurring  whenever  the  fuel -air  ratio  in  the 
colder  fan  duct  airstream  exceeds  some  critical  level.  This  critical  level  varies  with  the  exact 
geometry  of  the  flameholder  and  the  operating  conditions  of  the  augmentor. 

The  nature  of  the  feedback  mechanism  which  drives  the  instability  is  the  classical  pressure 
and  velocity  response  of  the  heat  release  mechanism  causing  variations  in  the  average  augmentor 
exit  temperature.  These  variations  result  in  oscillations  in  average  nozzle  inlet  pressure  through 
the  choked  flow  relationships  for  a constant  mass  flowrate.  This  results  in  further  oscillations  in 
the  flow  conditions  around  the  sprayring  and  flameholder  region.  The  relative  amplitudes  and 
phase  angles  are  functions  of  the  geometry  and  operating  conditions  of  the  augmentor  and  fan 
duct. 


This  program  was  directed  toward  the  analysis  and  modeling  of  the  heat  release  process  in 
terms  of  physical  geometry  and  actual  operating  conditions.  The  analysis  and  modeling  of  the 
response  of  the  augmentor  .system  to  the  sen.sitivity  of  the  heat  release  process  was  done  in  a 
companion  program  under  Contract  F3.'1615-7t-C-2024. 

The  analysis  and  model  for  the  heat  release  rate  were  .structured  to  rely  as  heavily  as 
possible  on  analyses  of  the  physical  processes  of  the  fuel  preparation,  flame  stabilization  and 
flame  propagation.  In  this  manner  it  was  felt  that  maximum  utility  of  the  model  would  be 
realized.  All  too  frequently  combustion  stability  models  are  generated  with  built-in  correlation 
constants  which  are  generic  to  one  form  of  combustor  .miy.  Whenever  a geometry  variation  is 
performed  there  is  no  guideline  for  the  required  change  in  those  constants,  and  the  utility  of  the 
model  is  limited  until  sufficient  experience  is  obtained.  In  this  model  format  we  have  attempted 
to  remove  that  restriction. 

The  mechanism  of  response  of  the  heat  release  process  to  variations  in  pressure  and  velocity 
over  the  sprayring  and  flameholder  region  is  described  bylow  with  reference  to  Figure  2. 

• A variation  in  pressure  and  velocity  over  the  fuel  injection  spraybar  results  in 
a variation  in  the  average  air  mass  flowrate  and  thus  local  fuel-air  ratio,  since 
the  fuel  flowrate  is  essentially  constant. 

• At  the  flameholder,  the  pressure  and  velocity  oscillation  results  in  a variation 
in  the  rate  of  formrtion  of  the  flameholder  wake  fuel-air  ratio  from  the 
available  liquid  phase  fuel. 

• The  pressure/velocity  oscillation  also  results  in  a variation  in  the  level  of  the 
wake  reaction  efficiency  at  the  above  level  of  fuel-air  ratio.  This  results  in  a 
variation  in  the  rate  of  initial  flame  propagation  into  the  approach  stream 
fuel-air  mixture. 
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• The  pressure/ve!(x-ity  oscillation  also  results  in  a variation  in  the  rate  of 
transverse  flame  spreadiOfC  at  the  local  flame  speed  due  to  the  influences  of 
local  pressure  and  turbulence  on  flame  speed. 

• The  final  result  is  a variation  in  the  average  value  of  the  nozzle  inlet 
temperature  due  to  the  combined  effects  of  the  previous  responses. 

This  variation  occurs  at  a time  delay  equal  to  the  transport  delay  between  the  sprayrings  and  the 
nozzle. 

This  proce.ss  includes  two  major  mechanisms  which  cause  a response  in  nozzle  inlet 
temperature  lO  oscillations  in  augmentor  inlet  pres.sure  or  velocity.  The  variation  in  local  fuel-air 
ratio  (i.e.  axial  variation  at  fixed  overall  level)  causes  an  axial  variation  in  heat  release  which  is 
felt  as  a temporal  oscillation  in  nt  zle  inlet  temperature.  Imposed  onto  this  variation  is  the 
variation  in  axially  local  heat  reler  e.  at  the  value  of  fuel-air  ratio,  caused  by  response  of  the 
flame-holding  process  to  the  pressure  or  velocity  oscillations. 

The  combustion  model  required  must  thus  he  able  to  relate  the  ultimate  average  of 
augmentor  exit  temperature  to  the  local  values  of  operating  conditions  at  the  flameholder  and 
fuel  injection  source  as  well  as  to  the  geometry  of  the  augmentor  and  flameholder. 

b.  Combustion  Model  Philosof  y 

“...  it  is  only  rational  to  be  in  with  an  idealization,  a simplification,  even  though 
we  know  that  to  be  ‘wrong’.  We  begin  with  that  as  a first  approximation,  then 
correct  for  major  discrepancies,  then  for  minor  di.screpancies,  then  for  very  minor 
discrepancies,  and  so  on.  Little  by  little  we  approach  a (possibly  unattainable) 
real  ‘truth’  and  in  the  prc^ess  develop  a precision  as  tight  as  neces.sary  for  our 
purposes.” 


Isaac  Asimov,  1969 

The  basic  framework  of  the  model  will  be  described  with  references  to  Figure  2.  This 
illustration  may  be  visualized  as  representing  a two  dimensional  sector  of  a flameholder  array.  A 
full  augmentor  analysis  would  require  a multistreamtube  solution  to  represent  the  full 
flameholder. 

The  augmentor  inlet  conditions  are  known  in  terms  of  a”ea,  pressure,  velocity,  temperature, 
turbulence  intensity  and  vitiation  level  (if  any).  The  required  overall  fuel-air  ratio  and  thus  fuel 
flowrate  from  the  spraybar  is  .specified  in  terms  of  flowrate  versus  pressure  drop  and  droplet  size 
and  distribution  also  versus  pressure  drop.  Once  specified,  the  flowrate  of  fuel  or  overall  fuel-air 
ratio  and  duct  static  pressure  will  define  a droplet  s’ze  distribution. 

At  this  point  of  injection  a simple  enthalpy  balance  is  performed  to  evaluate  the  percentage 
of  the  fuel  which  vaporizes  due  to  the  adiabatic  throttling  process  of  injection  from  the  high- 
pressure  sprayring  into  the  low-pressure  augmentor.  For  this,  the  fuel  properties  and  fuel 
temperature  in  the  sprayring  must  be  known.  The  fuel  flowrate  which  remains  in  the  liquid  phase 
is  placed  into  5 or  10  fuel  droplet  size  groups  which  represent  equal  mass  flowrate  distributions 
of  the  spray  distribution  curve  for  the  particular  sprayring  injector. 
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These  fuel  droplets  are  a’lowed  to  accelerate  with  the  airstream  tou^'rrls  'he  nanichi-hlt  r 
with  concurrent  droplet  vaporization.  The  vaporization  analysi.-  utilizes  a inn cd  < (uived inn 
model  for  combined  heat  and  mass  transfer.  A finite  difference  scheme  intenraies  the  accelerat  inu 
vaporizing  droplet  lifetime  until  either  the  llameholder  plane  is  reached  or  the  droplet  diameter 
goes  to  zero.  This  is  performed  for  each  droplet  size  group.  The  droplets  are  treated  as  sph.erical 
for  the  acceleration  and  vaporization  analysis. 

It  is  at  this  location,  the  flameholder  leading  edge,  that  the  comhiistion  model  riev dates  from 
classical  flameholder  or  combustion  analyses.  1'he  classical  e.cplanation  for  the  mec  hanism  hv 
which  a flameholder  functions  is  that  the  wake  of  the  bluff  body  serves  as  a volume  tdr  the 
reaction  of  the  vapor  phase  mainstream  fiicl-air  mi.vt  re.  Ignition  of  the  high-spi'ed  mainstream 
mixture  is  achieved  by  the  hot  wake  reaction  products  in  the  shear  layers  aft  of  the  llameholder. 
The  degree  of  reaction  achieved  in  the  wake  volume  is  dependent  on  the  ratio  of  a\  erage  residence 
tin.e  to  average  reaction  rate.  The  reaction  rate  depends  on  the  entering  vapor  fuel-air  ratio  and 
operating  pressure  and  temperature,  i.e.  kinetic  rate  constants. 

The  problem  is  that  in  a relatively  cool  airstream,  such  ns  the  fnn  duct  e.vit 

temperature,  and  overall  fuel -air  ratios  which  typically  react  wcdl  the  dc’grcc  of  circ;plet 
vaporization  is  very  low.  .So  low.  in  fact,  that  the  vapor  phase  fuel-air  ratio  at  the-  tlnn.eholdcr  is 
well  below  the  lean  flammability  limit  for.iP-type  fuels.  Howevc'r.  these  conditions  do  produce 
stable  flames  in  actual  burners.  Obviously,  there  is  some  addiiional  mechanism  by  which  vapor 
phase  fuel  is  generated  and  mixed  into  the  flameholder  wake. 

The  mechanism  which  has  been  identified  for  this  process  is  tlu'  formation  cif  a licpiid  film 
on  the  surface  of  the  flameholder  by  impingement  of  the  fuel  dro[)!e's.  This  surface  film  is 
partially  vaporized  hy  the  heat  flux  from  the  hot  wake  through  the  surlac  <■  of  the'  flame' holder  into 
the  film.  'I'his  vaporized  fuel  ent^'rs  the  recirculation  zone  through  the  .shear  layer  and  providc-s 
the  bulk  of  the  vapor  fuel  for  the  wake  reaction  process. 

For  the  purposes  of  the  model,  the  process  of  fuel  collection  on  the-  surfac-e  of  the  I'lamc-holdcr 
is  analyzed  as  spherical  droplet  trajectories  through  the  flow  field  a.^  it  moves  around  the 
flameholder.  The  analysis  is  performed  for  each  of  the  initial  droj)  size  gnnijis  utilizing  the 
diameter  which  remains  after  the  droplet  has  experienced  partial  vnpi'r’zntion  hc-tveen  the 
sprayhar  and  the  flameholder.  This  analysis  defines  the  fuel  collection  r.Je  onui  the  surface  of  I he 
flameholder,  when  integrated  over  the  droplet  size  groups. 

Once  the  film  is  established  on  Ihe  surface,  the  rate  of  vaporization  is  analszed  by  assuming 
a Nusselt  number  form  of  forced  mass  transfer  driven  by  the  vapor  jiressure  of  the  film.  A 
combined  film  heating  and  vaporizatior.  solution  is  performed  as  a finite  dif  ference  analysis  along 
the  surface.  The  heat  flux  is  evaluated  from  a wake  film  coefficient  and  ;m  assumed  wake 
temperature.  Since  the  wake  temperature  is  a function  of  the  fuel  vaporization  rate  and  the  wake 
reaction  level,  an  iterative  solution  is  required  between  this  initial  guess  of  wake  temiierature  and 
the  calculated  value. 

As  the  finite  difference  solution  priKeeds,  a certain  amount  of  litjuid  fuel  accumulates  which 
has  experienced  heating  but  no  vaporization.  This  fuel  is  lost  to  the  near  field  stabilization 
process  but  does  enter  into  the  far  field  reaction  during  turbulent  flame  sijreading.  'I'he  energv 
required  to  rai.se  this  fuel  from  i*s  collection  temperature  to  the  final  film  temiierature  represents 
a heat  loss  from  the  wake  reaction. 

For  the  evaluation  of  the  wake  reaction,  ihe  wake  is  treated  as  a well-stirrcd  chemical 
teactor  operating  with  gaseous  fuel.  To  perform  this  calculation,  the  volume  and  mass  influx  rate 
must  be  known.  These  values  are  produced  from  publi.shed  data  on  the  relative  size  and 
recirculation  rates  behind  bluff  body  stabilizers  as  functions  of  geometry  and  operating 
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conditions.  (Specific  source  references  v/ill  be  given  as  the  analyses  are  developed  in  a later 
section  of  this  report.)  The.se  data  have  been  reduced  to  a series  of  empirical  functions  which 
relate  wake  volume  and  t irculation  rate  to  flameholder  blockage,  geometry,  approach  f.ow 
conditions,  and  turbulence  level. 

The  reaction  efficiency  in  this  known  sire  wake  reactor  is  analyzed  assuming  that  the 
process  proceeds  as  a single  step,  .second  order  reaction  with  the  inefficiencies  represented  as 
unreacted  CO  for  lean  operation  or  unreacted  fuel  for  rich  operation.  The  entering  mass  flowrate 
is  balanced  against  the  mass  consumption  rate  at  an  unknown  final  level  of  efficiency.  A straight- 
forward solution  for  the  reaction  efficiency  and  thus  wake  temperature  proceeds.  Any  external 
heat  gain  or  loss  nechan'sms  are  used  to  adjust  the  reaction  rate  and  increase  or  reduce  the  final 
reactor  efficiency. 

At  the  known  level  of  entering  air  and  droplet  vaporized  fuel  flowrate  plus  surface  vaporii.ed 
fuel,  the  iteration  is  performed  as  follows: 

Recirculation  rate  calculated 
Wake  temper.Hture  assumed 
Surface  vaporization  calculated 
W'ake  fuei-air  ratio  thus  known  , 

Wake  reaction  level  at  this  f/a  calculated 
Wake  temperature  thus  known 
Iterate  back  through  surface  vaporization 

At  this  point  the  wake  conditions  are  known,  and  the  fuel-air  ratio  distribc'ion  around  the 
flameholder  is  known.  The  turbulent  flame  spreading  into  the  free-stream  is  initiated  in  the  shear 
layer  by  the  hot  wake  products.  The  degree  of  perfect  initiation  depends  on  the  excess  thermal 
energy  available  from  these  pnxlucts.  Flame  initiation  is  a go  or  no-go  phenomenon  and 
statistical  in  nature,  since  a lack  of  perfect  initiation  physically  results  in  local  regions  along  the 
surface  of  zero  ignition.  As  the  temperature  of  the  wake  products  is  reduced,  the  percentage  of  the 
area  which  fails  to  ignite  increases.  As  these  regions  increase  in  number  and  size,  a greater  portion 
of  the  heat  liberated  in  those  areas  which  do  ignite  is  transferred  into  adjacent  unlit  regions  rather 
than  into  transverse  flame  propagation.  Since  the  model  is  based  on  uniform  flame  initiation,  this 
process  is  observed  a?  a slower  rate  of  flame  spreading.  To  produce  this  re.sult,  the  model  relates 
the  flame  speed  to  the  ideal  turbulent  value  and  the  wake  efficiency  level. 

The  form  of  the  flame  spreading  model  is  a thin  sheet  flame  front  propagating  into  a fuel- 
air  mixture  at  a velocity  which  was  accelerated  by  the  flameholder  blockage.  The  specific  analysis 
for  the  turbulent  flame  speed  follows  Karlovitz,  where  the  turbulent  flame  speed  is  related  to  the 
laminar  flame  speed,  at  the  approach  conditions  and  the  local  fuel-air  ratio,  and  the  local  value 
of  turbulent  velocity.  This  latter  term  is  evaluated  ■ fijm  the  turbulence  generated  by  the 
flameholder  and  axially  decayed  in  a 10  L/D  length,  based  on  the  width  of  the  effective  jet  flow 
between  the  flameholders,  as  a function  proportional  to  x-'  i.  The  final  value  after  decay  is  the 
free-stream  intensity. 

An  additional  term  is  added  to  the  calculated  value  of  the  turbulent  flame  speed  to  account 
for  the  sensitivity  of  the  reaction  rate  to  integrated  efficiency  and  flame  self-turbulence.  These 
terms  are  modeled  as  a multiplier  whose  value  depends  on  the  local  efficiency  and  has  a value  of 
1 at  zero  and  lOOfr  efficiencies.  The  peak  value  is  2 at  5Cfr  efficiency. 

The  transverse  position  of  the  flame  front  is  found  by  a finite  difference  integration  of  the 
local  flame  front  into  the  approach  flow.  Due  to  the  sustaining  effect  of  the  local  heat  release,  the 
approach  velocity  is  axially  retained  at  the  value  accelerated  by  the  flameholder. 


\ 


12 


Tlie  overall  efficiency  is  found  at  the  nozzle  entrance  through  integration  of  the  flame 
spreading  heat  release  evaluated  from  the  preceding  analyses. 

c.  Combustion  Model  Program 

The  combustion  analyses  are  computerized  into  one  cohesive  program  which  combines  the 
models  for  the  various  processes  discussed  previously  under  Combustion  Model  Philosophy.  The 
computerzation  approach  selected  is  one  which  utilizes  separate  subroutines  for  each  of  the 
process  sections  of  the  model.  Although  this  approach  consumes  slightly  more  computer  run  time 
than  a full  step-by-rtep  calculation  procedure  it  was  .selected  for  two  reasons; 

1.  It  allows  easier  alterations  to  each  separate  process  should  changes  in  any  of 
the  local  analyses  be  desired,  e.g..  a revi.sed  reaction  kinetics  mf)del.  This 
reduces  the  chances  for  errors  due  to  unforeseen  interactions  between 
sections  of  the  program. 

2.  The  it.  ve  steps  required  between  the  wake  analysis  and  the  surface 
vaporiza!  n analysis  is  facilitated  if  they  are  isolated  in  the  program. 

Written  in  this  manner,  a single  inlet  condition  analysis  requires  only  about  10  seconds  of 
CPU  execution  on  the  IBM  .“ITO  168  computer.  The  required  input  consists  of  a description  of  the 
geometry  of  the  streamtube: 

Bloc’  age  ratio 
Flameholder  width 
Flameholder  apex  angle 
Spraybar  to  flameholder  spacing 
Flameholder  to  exhoast  spacing. 

Those  parameters  which  describe  the  operating  conditions  are: 

Inlet  pressure 
Inlet  velocity 
Inlet  temperature 
Inlet  turbulence  intensity 
Fuel-air  ratio 
Spraybar  fuel  pressure 
Spraybar  fuel  temperature 
Fuel  type 

Vitiation  fuel-air  ratio. 

These  inputs  define  the  required  flamework  for  the  analysis.  The  pt'ogram  organization  is 
shown  in  Figure  3.  The  various  subrt/utines  perform  the  analyses  as  follows: 

INJECT  This  performs  the  spray  formation  and  throttling  process 
vaporization  analyses. 

ACCEL  This  evaluates  the  forced  vaporization  and  acceleration  of  the 

droplets  between  the  spraybar  and  flameholder.  Done  one.*  per 
droplet  size  group. 

COLLECT  This  evaluates  the  collective  rate  of  the  liquid  droplets  onto  the 
surface  of  the  flameholder.  Performed  once  per  droplet  size 
group. 
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RECIRC 

This  evaluates  the  size  of  the  recirculation  zone  wake  behind  the 
flameholder  and  the  rate  of  entry  of  gaseous  comjroiients  into 
this  volume. 

B.3 

This  evaluates  the  rate  of  vaporization  of  the  liquid  film  from  the 
surface  of  the  flameholder. 

WAKE 

For  the  results  of  RE(TRC  and  B.1.  this  solves  for  the  wake 
reaction  efficiency  and  temperature. 

FLAME 

This  performs  the  solution  for  the  transverse  penetration  of  the 
turbulent  flame  from  the  woke  shear  layer  into  the  free-stream 

mixture. 

Currently,  the  input  and  output  reflect  the  requirement.*  to  analyze  a sinple  case  streamline 
analysis.  The  output  defines  the  intermediate  results  of  the  various  subroutines  as  well  as  the 
heat  release  profile  between  the  flameholder  and  the  exhaust  nozzle.  The  program  is  capable  of 
multiple  case  execution  limited  only  by  run  time.  This  capability  would  he  utilized  for  multiple 
streamtube  analysis  and  parametric  studies.  At  pre.sent.  there  is  no  ^aphics  output  capability, 
but  it  could  be  added. 

d.  Analyses  and  Results 

The  combustion  model  as  formulated  for  the  rumble  mechanism  essentially  consists  of  two 
parts.  These  are  a compositional  analysi.s  and  a reaction  analysis.  The  compositional  analysis 
defines  the  manner  in  which  a combustible  fuel-air  mixture  is  generated  from  the  injected  liquid 
fuel.  The  reaction  analysis  defines  the  manner  in  which  this  mixture  actually  bums  to  produce 
a final  level  of  combustion  efficiencv. 

The  anal>sis  for  the  compositional  portion  is  outlined  below  for  the  basic  ca.se  of  a two- 
dimensional  duct  *•  ifh  a bluff  body  stabilizer  and  liquid  fuel  injection  source.  For  ease  of 
development  of  the  .solutions,  the  following  as.suniptions  were  u.sed: 

• The  injected  liquid  fuel  forms  a homogeneou.siy  di.spersed  spray  which  fills 
the  full  cross  section  of  the  duct. 

• The  approach  flow  field  is  uniform  in  the  transverse  direction,  i.e.  trunsverse 
uniformity  of  velocity,  pressure,  and  temperature. 

— Under  these  assumptions,  the  form  of  the  compositional  analysis  was  developed.  The 
schematic  «/  the  proces.ses  i'  shown  in  Figure  4. 

A certain  portion  of  the  injected  fuel  flowrate  experiences  flash  vaporization  during  the 
injection  proces.s  and  forms  a uniform  flowrate  of  vapor  phase  fuel.  This  fuel  flowrate  is; 

Wf,  =-  di>i.  (1) 

The  remaining  portion  of  fuel  goes  into  a liquid  droplet  spray  whch  undergoes  spray  vaporization. 
The  vapor  flowrate  produced  by  this  process  is: 

Wf,  =•  diE<l-dir)Wf.  (2) 
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Main  ) — I Input 


Figure  3.  Ducted  Flame  Model  Linkage  Map 


Heating  and  Droplet  Vaporization, 


Figure  4.  Schematic  of  Wake  Mixture  Formation 


Combining  these  processes  into  one  influence  coefficient  we  may  write: 

Wr,  = fitW,.  (3) 

We  continue  to  utilize  the  assumption  that  the  vaporized  fuel  is  homogeneously  spread  as  it  is 
produced.  The  next  processes  which  are  evaluated  are  those  concerned  with  the  formation  of  the 
liquid  film  vaporization  rate. 

At  the  flameholder  plane,  the  mixture  consists  of  a flowrate  of  air  (w„),  vapor  fuel  (wr^)  and 
liquid  fuel  (w,,)  defined  as: 

w,  = p»VA  (4) 


Wfv  = ^>WfT  = w,  (6) 

wr,  = w,.^  - w,^  = (1-/3,)  w,.  (7) 

This  leaves  (1-/3,)  liquid  fuel  percentage  available.  Flowing  through  the  flameholder  area, 
then  we  have: 

w.p  = Fw,  = F/j,  VA  (8) 

wr,r  = (1-/9.)  flw.F  (9) 

w,,,.  = /3,fl  w.F.  (10) 

The  liquid  and  vapor  fuel-air  ratios  at  this  point  are: 

, /9i9w,F 


(f/a),  = (1-/3,)  e.  (12) 

A portion  of  the  liquid  fuel  is  now  collected  and  vaporized  from  the  surface.  The  collection 
rate  of  liquid  fuel  is: 

Wf  ^ = /JiWpr  = /3,(l-/3.)  <F;  (13) 

of  this  amount,  ff,  percentage  vaporizes  from  the  surface  and  recirculates: 

w = /3,w,,_.  4=  fit  did-/?,)  e w.F.  (14) 

The  prevaporiz^d  fuel  recirculates  at  a rate  of  K,  percent  which  passes  the  blocked  area: 

w = Wf,,.  K,  - F K,.  (15) 

'l 

Thui  the  total  recirculated  vapor  fuel  flowrate  is: 

1 

Wf,,  = d,9w.FK,  -t-  (l-/3,)d,di9w.F.  (16) 
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The  air  recirculation  rate  is  K,  of  tlie  flowrate  through  the  blocked  area; 

W.,„  = w„rK,.  . (17) 

From  Equations  16  and  17,  the  wake  fuel-air  ratio  in  vapor  phase  is: 


<t>  = 

Wf 

'vw 

( 18) 

<p  = 

0,0  + {l~0,)O 

K, 

(19) 

A- 
e ~ 

/3.  + 

K, 

(20) 

Equation  20  relates  the  wake  vapor  phase  fuel-air  ratio  to  the  overall  fuel-air  ratio  and. 
through  the  influence  coefficients,  to  the  geometric  and  aerothermodynamic  situation  under 
scrutiny.  Each  of  these  may  now  be  evaluated  from  the  available  body  of  combustion  literature 
to  determine  the  overall  form  of  the  compositional  model. 

A comment  on  the  implication  of  Flquation  20  is  in  order  here.  This  equation  determines  the 
mixture  in  the  wake  which  then  undergoes  reaction  to  provide  the  ignition  source  for  the 
mainstream  reaction.  As  with  any  stirred  chemical  reactor,  the  mosc  efiicient  process  for  a fixed 
volume  occurs  near  the  stoichiometric  fuel-air  ratio.  Since  the  degree  of  tins  reaction  efficiency 
strongly  influences  the  overall  augmentor  efficiency,  the  fuel-air  ratio  should  be  kept  near 
stoichiometric. 

If  we  introduce  the  following  values  from  the  analyses  to  be  presented  later: 


= 0.20 

02  = 0.75 

02  = 0.50 

K.  = 0.25 


we  find  that  Equation  20  yields; 


= 0.20  -f-  (0.8) 
u 


(0.75)(0.,50) 

0.25 


A - ISO. 

Thus,  due  to  the  relative  fuel  concentrating  effect  of  the  flameholder,  the  wake  is  40''i  richer  than 
the  overall  system.  If  the  overall  fuel-air  ratio  is  0.0.50,  the  wake  is  over  stoichiometric  at  0.070. 
Any  further  increase  in  fuel  flowrate  results  in  a drastic  decrease  in  the  wake  reaction  efficiency 
due  to  the  increase  in  wake  f/a. 

The  second  portion  of  the  analyses  is  the  wake  reaction  efficiency  and  turbulent  flame 
spreading.  These  analyses  will  be  developed  in  the  following  sections. 
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(1)  Fuel  Vaporization  Before  Flameholder 


As  developed  in  the  preceding  analysis,  the  percentage  of  the  liquid  hydrocarbon  fuel  which 
vaporizes  prior  to  reaching  the  flameholder  consists  of  two  mechanisms: 

The  forniulion  of  vapor  due  to  the  throttling  process  of  injection  (diT-)- 

The  evaporation  of  fuel  droplets  due  to  forced  convection  (/Sib)- 

Thus  we  define  /3i  as: 


= /3.x  + (21) 

The  process  of  flash  vaporization  during  the  process  of  injection  is  treated  as  an  adiabatic 
flow  from  a region  of  high  pressure  and  moderate  temperature  (within  the  sprayring)  to  a low- 
pressure  area  (in  the  augmentor).  This  flow  situation  is  analogous  to  the  expansion  valve  process 
and  is  evaluated  from  knowledge  of  the  final  properties  and  operating  conditions. 

The  process  is  evaluated  assuming  adiabatic  expansion  from  known  levels  of  fuel  pressure 
and  temperature  within  the  sprayring  to  a kn.ywn  pressure  level  in  the  augmentor.  If  the  sprayring 
pressure  is  sufficient  to  maintain  the  liquid  fuel  as  a saturated  liquid,  the  enthalpy  level  is 
defined  for  the  given  fuel  type  as  a function  of  the  sprayring  fuel  temperature  only. 

For  adiabatic  injection,  the  final  mixture  enthalpy  equals  the  saturated  liquid  enthalpy  at 
T.  The  mixture  enthalpy  and  a known  static  pressure  will  define  the  quality  of  the  injected  fuel. 
This  process  is  readily  evaluated  from  the  enthalpy  diagram  of  the  particular  fuel  type.  Enthalpy 
diagrams  for  JP-4  and  JP-5  are  shown  in  Figures  5 and  6. 

An  expanded  diagram  of  the  adiabatic  process  for  JP-4  is  shown  in  Figure  7.  This  example 
is  for  injection  of  250“^  fuel  into  a 7.5  psia  environment.  The  parametric  results  for  JP-4  fuel  are 
shown  in  Figure  8. 

The  remaining  liquid  fuel  is  partitioned  into  the  droplet  size  groups  on  a 5 or  10  group  equal 
flowrate  basis.  For  this,  the  droplet  size  distribution  must  be  known.  Due  to  the  large  number  of 
spray  systems  available  to  the  augmentor  designer,  no  attempt  was  made  to  model  this  droplet 
formation  process.  Rather,  a distribution  function  is  built  into  the  program  which  describes  the 
size  distribution  as  a function  of  the  fuel  pressure  drop.  The  curve  used  in  the  program  represents 
the  droplet  distributions  for  variable  area  pintle  spraybars.  Figure  9 shows  this  function. 

If  5 size  groups  are  to  be  used,  each  one  represents  20%  of  the  liquid  flowrate.  The  sizes  used 
in  the  analysis  would  be  the  10,  30,  50,  70  and  90%  diameters.  Once  these  sizes  are  known,  the 
program  performs  the  solution  to  the  forced  convection  droplet  vaporization  between  the 
spraybar  and  the  flameholder. 

The  basis  for  the  droplet  vaporization  solution  is  the  form  of  the  Nusselt  number  devised  foi 
spherical  droplets  by  Ranz  and  Marshall  (Reference  1)  and  subsequently  improved  by  Rreim  and 
Heidmann  (Reference  2).  The  assumed  forms  are: 

Nuh  = 2 + 0.6  Re*'*  Pr*'*  (22) 

Num  = 2 + 0.6  Re*'*  Sc*'* 
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(23) 


for  heat  and  mass  transfer.  For  the  assumption  of  Lewis  No.  = 1,  i.e.,  Pr  = Sc,  these  become 
identical.  The  fuel  droplet  mass  efflux  rate  is  calculated  from: 


(24) 


Nu  Dv  MW 

RdfT.  • 


(25) 


During  the  transient  droplet  heating  period  between  injection  at  T = Tf^  and 
achievement  of  the  wet-bulb  temperature,  the  droplet  temperature  is  evaluated  from  the  heat 
input  and  the  diffusion  driven  mass  efflux.  The  difference  between  the  heat  requirea  ‘o  generate 
the  mass  efflux  and  the  heat  which  actually  reaches  the  surface  of  the  droplet  is  assumed  to  go 
towards  alteration  of  the  droplet  bulk  temperature.  This  assumption  is  essentially  a statement 
that  the  droplet  internal  circulation  is  sufficiently  rapid  compared  to  the  thermal  input  that 
significant  droplet  radial  thermal  gradients  do  not  exist.  Although  this  assumption  has  been 
questioned,  most  noteably  by  Strahle  (Reference  3),  any  errors  introduced  are  quite  secondary  to 
errors  in  evaluation  of  the  droplet  film  transport  properties.  Sufficient  accuracy  and  significant 
numerical  simplification  is  realized  by  this  method. 

In  evaluation  of  the  net  heat  flux  into  the  liquid  surface,  the  mass  efflux  blocking  term 
introduced  by  Preim  (Reference  2)  is  used.  This  term  evaluates  the  loss  in  net  flux  due  to  the 
vapor  superheating  which  occurs  as  the  evolved  fuel  vapor  achieves  thermal  equilibrium  at  the 
free-stream  gas  temperature.  The  transient  solution  proceeds  as  follows,  after  evaluation  of  the 
fuel  mass  efflux  rate,  w,  from  Equation  24. 


The  thermal  film  coefficient,  hf,  is: 
t k Nu„ 


(26) 


where  Nuh  is  from  Equation  22. 

The  net  heat  flux  to  the  droplet  liquid  surface,  q,  is  calculated  with  allowance  for  the 
thermal  blocking  due  to  vaporizing  fuel  heating  (/3): 

q = h,A.(T,-T, ) p (27) 

where 

^ (28) 

and 


z = Cpv  w/x-kd,  Nuh.  (29' 

The  net  amount  of  heat  which  is  available  for  sensible  heating  of  the  droplet  is  the  net  liquid 
surface  flux  minus  the  latent  heat  required  to  generate  the  vapor  mass  flux: 

Aq  = q - WvA.  (30) 
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Fuel  Injection  Temp  - °F 

FI)  1141144 

Figure  7.  Constant  Enthalpy  Fuel  Injection  of  J P-4  Fuel 


Figure  8.  Throttling  FYocess,  % Vaporized  vs  Fuel  Temperature  (JP-4  Fuel) 


Figure  9.  Typical  Spray  Distribution 


This  net  flux  is  used  to  raise  the  liquid  bulk  temperature,  T,  as; 


dT,  _ .Iq  _ Aq 

dt  m.  Cp,  4/3  » {d,/2)*  p(  Cp,  ' 


(31) 


This  change  in  T may  be  positive  or  negative  depending  on  the  relative  values  of  w and  q. 

The  wet-bulb  temperature  is  defined  as  that  liquid  droplet  temperature  where  the  thermal 
net  input  is  just  sufficient  to  generate  the  mass  ffflux  at  that  temperature  and  vapor  pressure.  As 
such,  it  is  a function  of  effective  Nusselt  number  and  local  static  pressure  and  temperature. 
Coincident  with  this  transient  heating  and  vaporization,  the  liquid  droplet  is  being  accelerated  by 
the  faster  free  stream  gas  velocity.  The  acceleration  is  evaluated  from  the  standard  equation: 


dV, 

dt 


3 Cd  Pa 

4 d,  p. 


(V.-V,)* 


(32) 


Since  the  Reynolds  number  used  in  the  Nusselt  number  formulation  is  defined  as; 


Re  = 


Pad,  (V.-V,) 


(33) 


a simultaneous  solution  of  Equations  31  to  33  is  required.  In  the  computer  analysis,  this  is 
accomplished  by  a finite  difference  solution  utilizing  small  time  intervals.  The  ordinary 
differential  equations  are  rewritten  as  delta  terms,  e.g.. 


dt 
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becomes 


M • 

For  each  time  increment,  the  initial  values  are  used  to  calculate  Re.  Nuh.  N’um.  pv.  etc. 
Equations  24  to  .31  are  solved  which  yields  Awv  and  AT  . Equation  32  is  solved  to  yield  AV’  . At 
the  end  of  this  step,  w»  is  incremented  by  AWv.  T by  AT  and  V by  AVMTie  average  value 
of  V over  this  At  yields  a delta  axial  travel  distance. 

The  procedure  is  repeated  in  small  time  steps  until  either  the  total  axial  trawl  exceeds  the 
spraybar-to-flameholder  separation  distance  or  the  liquid  is  fully  vaporized.  At  the  start  of  each 
new  time  increment,  the  input  values  reset  to  reflect  the  effect  of  the  previous  step;  e.g..  pv  is 
evaluated  at  T<  + AT  , d,  is  evaluated  to  reflect  less  liquid  mass,  etc. 

Some  typical  results  for  the  vaporization  portion  of  the  analysis  are  presented  in  Figures  10 
to  12.  For  these  results,  the  droplet  acceleration  portion  of  the  analysis  was  deleted  to  allow 
solution  on  a limited  size,  real-time  access  computer  terminal.  Thus,  the  relative  veliH-ity  is  held 
constant  and  the  results  represent  the  upper  limit  on  transient  heating  and  mass  e jx.  Figure 
10  shows  the  transient  bulk  liquid  temperature  for  a 100  micron  droplet  in  a 400°F.  250  ft/sec  gas 
stream.  The  decrease  in  temperature  from  the  initial  SO^F  injection  temperature  for  the  0.5 
atmosphere  pressure  reflects  the  fact  that  the  initial  mass  efflux  by  convection  c reeds  the 
available  net  heat  input.  The  wet-bulb  temperature  as  a function  of  gas  static  pressure  and 
temperature  is  shown  in  Figure  11  for  the  same  100  micron  JP-4  droplet. 

Corresponding  with  the  wet -bulb  temperature,  the  equilibrium  ma.ss  efflux  for  a 1(X)  micron 
droplet  is  shown  in  Figure  12.  The  effect  of  pressure  and  temperature  is  clearly  shown. 


Figure  10.  JP-4  Droplet  Transient  Temperature 
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Figure  II.  JP-4  Fuel  Wet  Bulb  Temperature  vs  Pressure  and  Gas  Temperature 
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Figure  12.  JP-4  Forced  Convection  Vaporization 
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For  a situation  described  as  follows: 


V,  = 250  fps 

T.  = 275®F 

T = SO'-F 

P,  = 10  psia 

d,  = 25  to  125  microns. 

The  amount  of  JP-4  which  will  vaporize  in  a typical  di.siance  of  6 inc  hes  is  shown  in  Figure 
13.  These  results  are  typical  for  fan  duct  conditions  at  high  altitude  and  low  flight  Mach  number. 


Figure  13.  Vaporization  of  JP-4  Fuel  Droplets  I’s  Axial  Length 


(2)  FItmmholdpr  Fuel  Colltctlon  Efficioney 


The  output  from  the  ACCEL  subroutine  defines  a spray  of  droplet  diameters  at  a plane  wit  h 
the  flameholder  leading  edge.  The  diameters  are  defined  from  the  final  step  in  the  forced 
vaporizatk/n  analysis  as  is  the  droplet  axial  velocity.  The  spray  field  is  still  assumed  to  he 
homogeneous  in  the  transverse  direction  in  terms  of  size  and  volume  flowrate  distribution. 
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A*  the  flow  stream  approaches  the  stabilizer,  the  suspended  droplets  are  unable  to  fully 
follow  the  divergence  of  the  flow  strebmiines  and  a portion  of  them  impinge  on  the  bluff  body.  The 
evaluation  of  this  “capture  rate"  is  done  in  nondimensionalized  terms  as  fit'. 

Wic  = Wn  rd»  (34) 

(35) 


fuel  liquid  capture  rate 
fuel  liquid  flowrate  at  F/H  plane 
F/H  blockage 

dt  w Also  the  ratio  of  the  separation  of  the  droplet  capture  limit  streamlines  tc  the  F/H 
width; 


where 


rwfi 


Wfr  = 

w,  = 

r = 


w„  = P,V6*  (36) 

Wf.  * p,  V,  A 


Substituting  yields: 

A “ -fj-  (37) 

where  i*  - droplet  capture  width.  This  is  shown  schematically  in  Figure  14. 


Plane  A 


Figure  14.  Droplet  Capture  Schematic 
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This  annlysis  will  be  performed  once  for  each  sire  Krf’np  cxitinp  the  vnporirnf  ion  analysis. 
If  10  groups  are  used  initially,  each  sixe  represents  HK*  of  the  liquid  fuel  at  plane  A.  since  the 
smaller  sires  will  have  higher  vaporization  rates  between  the  injection  plane  and  plane  A.  The 
overall  collection  rate  is  thus: 


1 

m,  , A 


L 


ttti  i Adii , 


(W 


Evaluation  of  (i„  requires,  therefore,  evaluation  of  4*  for  any  value  of  droplet  diameter. 
Additionally,  the  streamline  divergence  rate  will  be  a function  of  the  flameholder  geometry  (aiiex 
angle,  blockage).  The  droplet  acceler  ition  rate  in  the  x and  y directions  will  be  a function  of  ihe 
flow  conditions  (pressure,  temperature,  velocity)  and  initial  fuel  droplet  velocity  (output  i.om 
0„  subroutine). 


An  earlier  analysis  for  collection  into  a cylindrical  bar  was  performed  by  Langmuir 
(Reference  4).  The  collection  efficiency  was  found  to  be  dominated  by  air  velocity  and  droplet 
diameter,  'fhe  influence  of  diameter  is  exceptionally  strong.  Typical  results  are  shown  for  vee 
gutter  flameholders  in  Table  1.  Going  from  25  to  75  micron  droplets  doubles  the  collection  rate 
from  ~30'f  to  ~70'r.  A typical  streamline  result  is  shown  in  Figure  15. 

The  effect  reinforces  the  dependence  of  the  low  air  temperature  stabilization  process  on  the 
surface  film  analysis.  As  the  air  temperature  is  reduced,  the  level  of  droplet  vaporization  is  also 
reduced  and  the  droplet  diameter  at  the  flameholder  is  increased.  This  not  onlv  increases  the 
liquid  flowrate  available  for  collection  but  also  increases  the  percentage  collected.  This  doubling 
effect  produces  a much  greater  dependence  of  the  wake  fuel-air  ratio  on  the  surface  film. 

TABLE  1 

COLLECTION  EPHCIENCY  VS 
DROPLET  DIAMETER 

Gutter  Apex  Blockage  Droplet  Collection 
Angle  Ratio  Diameter  Efficiency 
(deg^ (%)  (microns)  (%) 
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0.70 
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0.76 
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0.65 
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0.73 
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0.78 
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0.30 

40 
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60 

0.55 

40 

30 

76 

0.63 

40 

30 

100 

0.71 
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(3)  Fuel  Film  Surface  Vaporization  E¥aluatlon 


The  liquid  droplets  which  contact  the  surface  of  the  flameholder  are  assumed  to  generate  a 
uniform  liquid  film.  This  liquid  film  experiences  heat  transfer  from  the  hot  wake  through  the 
flameholder.  This  heat  is  partitioned  into  forced  mass  transfer  (latent  heat)  and  temperature  rise 
(sensible  heat).  See  Figures  16  and  17. 

The  liquid  film  vaporization  rate  is  used  to  calculate  the  surface  vaporization  efficiency, 
dt.  This  efficiency  is  the  percentage  of  the  collected  liquid  fuel  which  is  vaporized  from  the 
flameholder  surface  and  recirculated  into  the  wake  reaction. 

The  controlling  parameter  becomes  a forced  diffusion  process.  The  controlling  equation  may 
thus  be  written  as: 


Wfc  = Wv  + W,  + • 


Wv  = fcn(Num,  p„  A.) 
p,  = fcn(T, ) 

T,  s T,  c + Vi  AT 
q = w,X  + w,r  Cp  ATf 
q = fen  (Nuh,  T»,  T,,). 
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Figure  16.  Flameholder  Surface  Vaporization  Schematic 
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Figure  17.  Expanded  View  of  Region  A Heat  Transfer  Mechanisms 

In  these  equations,  the  collected  mass  flowrate,  is  thf*  result  of  the  B2  subroutine.  VVe  assume 
that  the  liquid  film  total  mass  does  not  change  with  time  unless  the  entry  flowrate  alters,  i.e. 


The  various  film  properties  are  functions  of  the  fuel  type  and  temperature.  Suitable  curves 
are  used  in  the  program  for  JP-4  and  JP-5  fuel. 

This  situation  is  a forced  mass  transfer  from  the  surface  represented  by  a Nusselt  number 
function  driven  by  the  film  vapor  pressure  and  a Nusselt  number  form  of  heat  flux  from  the  wake 
into  the  film. 

The  formulation  of  the  rate  of  liquid  fuel  vaporization  from  the  surface  of  the  flameholdcr 
has  been  programmed  as  a finite  element  solution  to  the  nonuniform  forced  diffusion  process.  The 
diffusion  process  is  evaluated  from  the  surface  Nusselt  number  for  mass  transfer: 


w,  = C.A.P.ln(^;-^) 


(45) 


C.  = 


Nu„  Dv  MW 
R AxT. 


(46) 


The  Nusselt  number  is  functionally  identical  for  heat  flux  and  mass  flux  when  the  Prandtl 
number  in  the  thermal  Nusselt  number  is  replaced  with  the  Schmidt  number  for  mass  transfer. 
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The  heat  flux  from  the  recirculation  zone  is  evaluated  directly  from  a Nusselt  num’jer 
formulation  for  recirculation  wakes  behind  bluff  bodies  as: 


h N 

Nu*  = — = 0.99  Re"  * P,"  ” 
km 


(47) 


For  a known  value  of  heat  flux  through  the  flameholder,  the  rate  of  surface  vaporization  is 
evaluated  by  breaking  the  liquid  film  into  10  subgroups.  In  each  subgroup,  the  vaporization 
process  is  assumed  to  be  controlled  by  the  vapor  pressure  corresponding  to  the  mean  liquid  film 
temperature  within  that  subgroup  (See  Figure  18). 


The  heat  flux  into  each  subgroup  is  responsible  for  two  physical  effects.  The  heat  is  used  to 
provide  the  latent  heat  of  vaporization  and  provide  sensible  heat  to  elevate  the  bulk  liquid 
temperature. 

For  the  initial  subgroup,  the  entrance  mass  flowrate  is  the  prorated  portion  of  the  fuel 
collection  rate.  Thereafter,  the  entry  flowrate  consists  of  the  prorated  collection  rate  plus  the 
unvaporized  liquid  fuel  from  the  preceding  subgroup. 


The  entrance  fuel  temperature  for  the  initial  subgroup  is  the  droplet  bulk  temperature  as 
evaluated  in  the  spray  vaporization  subroutine.  For  succeeding  subgroups,  the  initial  fuel 
temperature  is  evaluated  from  the  mixture  of  captured  fuel  at  droplet  temperature  and  liquid  fuel 
film  at  elevated  temperature. 

The  solution  in  each  subgroup  for  the  liquid  temperature  rise  and  vaporization  rate  requires 
an  iteration  process.  The  initial  guess  on  vaporization  is  evaluated  at  the  fuel  entry  temperature. 
The  latent  heat  required  to  accommodate  this  vaporization  is  subtracted  from  the  heat  flux  and 
the  excess  used  to  increase  the  bulk  film  temperature.  A new  vaporization  rate  is  calculated  at  a 
mean  fuel  temperature  between  entry  and  exit.  The  solutions  for  vaporization  and  temperature 
rise  are  repeated  until  convergence  occurs. 

This  process  of  finite  element  solution,  if  the  subdivisions  are  fine  enough,  allows  a 
relatively  simple  evaluation  of  the  vaporization  rate  for  nonadiabatic  walls  and  a nonuniform 
liquid  film  temperature.  It  may  be  made  as  exact  as  desired  by  decreasing  the  subgroup  size.  This 
may  be  very  desirable  for  solutions  to  multicomponent  fuels  (such  as  JP-4)  where  the  latent  heat 
of  vaporization  and  surface  vapor  pressure  are  strong  functions  of  the  bulk  liquid  temperature. 


The  partitioning  of  the  heat  flux.  Equation  43,  and  the  forced  diffusion  process.  Equations 
45  and  46,  may  be  combined  to  yield: 


q = w,^  (T,  ,-T, ,)  + X ( 


NUn,  Dy  MW 
RA.T. 


) A.p,ln  ( 


-Pi-) 
Pi-P»  / 


(48) 


where  Cp,  X and  p,  are  all  functions  of  Tf. 

The  functional  unknown  in  this  group  is  X t the  average  film  temperature.  Once  T is 
known,  the  mass  efflux  is  known  from  Equation  45.  The  form  of  the  efflux  response  to  T,  is  such 
that  a critical  film  temperature  exists,  when  pv  = p^  where  infinite  flux  exists.  For  JP-4  and  10 
psia  static  pressure  this  is  180®F.  A very  careful  analysis  is  required  to  iterate  to  a successful 
answer.  A typical  result  is  shown  in  Figure  19  for  JP-4  fuel  and  flow  conditions  of  250  fps,  275®F, 
10  psia. 
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Figure  18,  Finite  Difference  Solution  Procedure 
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Figure  19.  Surface  Vaporization  vs  Fuel  Collection  Rate  for  Two  Levels  of  Wake  Temperature 

(4)  FlamBhotder  Wake  Recirculation  Rate 

In  a manner  analogous  to  the  deHnition  of  the  collection  efficiency,  we  will  define  a 
recirculation  efficiency,  Ki.  This  efficiency  relates  the  mass  flowrate  into  the  wake  of  the 
flameholder  to  the  mass  which  flows  through  the  projected  blocked  area  of  the  flameholder: 

w,,  = p.V.NK,.  (49) 


This  recirculation  flowrate  is  the  mass  which  is  transferred  by  turbulent  diffusion  across  the 
free  shear  layers  aft  of  the  flameholder.  These  layers  form  the  boundaries  of  the  backflow  wake 
aft  of  the  bluff  body.  If  we  look  at  this  wake  volume  as  a homogeneous  region  with  mass  transfer 
across  the  boundaries,  the  recirculation  rate  may  also  be  written  as: 


w.,  = 


P.V„ 


(50) 


If  the  wake  volume  and  residence  time  may  be  evaluated  as  a function  of  the  geometric  and 
flow  variables,  then  the  recirculation  rate  may  be  evaluated  from  the  known  variables  without 
resorting  to  the  much  more  difficult  solution  to  the  effective  transport  across  the  shear  layers. 
This  approach  has  been  taken  by  the  majority  of  investigators  for  wake  processes  (References  5 
through  11). 


The  approach  relates  the  wake  volume  and  residence  time  to  the  controlling  aerodynamic 
-and  geometric  variables,  such  as; 


Blockage  ratio 

Vee  gutter  apex  angle 

Flow  Mach  number 

Pressure,  velocity,  temperature. 
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Once  these  are  known,  the  recirculation  rate  calculation  proceeds  as  follows: 

w,,  = p,  Vo/t  (51) 

Vo  = Cv(UD)(B/D)N*.  (52) 

If  we  nondimensionalize  the  residence  time  with  respect  to  velocity  and  characteristic  dimension, 
flameholder  width,  we  have: 


zYi 

N 


(53) 


V.  P.V, 
N t'  • 


Now,  from  Equations  52  and  54,  we  have: 

. _ p.V.Cv(L/D)(B/D)N 


(54) 


(55) 


From  this  and  Equation  49,  we  may  write: 

K,  = Cv(L/D)(B/D)(t')-«.  (56) 

The  value  of  the  shape  coefficient,  Cy,  was  determined  by  assuming  that  the  wake  was 
similar  to  a two-dimensional  ellipse  which  is  truncated  by  the  flameholder  at  its  forward  edge. 
The  value  used  is  0.80. 

The  data  presented  in  References  o through  11  were  reduced  to  a series  of  curves  which 
describe  the  effect  of  vanous  parameters  on  recirculation  rate.  The  functional  groupings  and  the 
corresponding  figures  are: 


L/D 

vs  }''lockage 

(Figure  20) 

L/D 

vs  Apex  Angle 

(Figure  21) 

L/D 

vs  Flow  Mach  Number 

(Figure  22) 

B/D 

vs  Blockage 

(Figure  23) 

B/D 

vs  Apex  Angle 

(Figure  24) 

t' 

vs  Blockage 

(Figure  25) 

r' 

vs  Air  Temperature 

(Figure  26) 

The  only  one  of  these  which  dc'^  not  follow  directly  from  the  cited  references  is  the 
temperature  dependence  of  the  recirculation  zone  residence  time.  The  data  available  were  taken 
over  an  extremely  limited  inlet  temperature  range  so  that  no  dependence  could  be  determined. 

The  shape  of  this  curve  was  drawn  by  assuming  that  the  forcing  function  for  the 
recirculation  mass  flowrate  is  the  turbulence  generated  at  the  recirculation  zone  boundary.  Since 
the  viscosity  of  the  free-stream  inhibits  this  formation,  the  curve  is  drawn  to  follow  the  increase 
in  viscosity  with  temperature  exhibitr  1 by  air. 
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L/D  (a)/L/D  (180)  Wake  L/D 
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Figure  20.  Wake  L/D  vs  Blockage  Ratio 
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Figure  24.  Effect  of  Apex  Angle  on  Wake  BID 
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Figure  26.  Influence  of  T,  on  Residence  Time 

For  a typical  analysis: 

V.  = 260  fps 
T,  = 275‘’F 
r - 0.25 

a = 60“ 

Mach  No.  * 0.189 
we  obtain: 

L/D  = 7.87 
B/D  = 1.55 
r'  = 36.9 

or: 

K,  = 0.265. 

Expre88e<f  in  direct  values,  at  10  psia  static  pressure  and  with  a 2-in.  flameholder: 

Wake  volume  = 39  in.*  per  inch  length 
Recirculation  rate  = O.C^  Ibm/sec  per  inch  length. 

The  total  duct  flowrate  for  this  case  is  0.51  tbm/sec  per  inch  length.  For  this  case,  6.7*7  of  the  free- 
stream  flowrate  is  recirculated  to  provide  ignition  for  the  remaining  93.397.  The  average  residence 
time  is  24.6  milliseconds. 

Before  going  on  to  the  reaction  kinetics  in  the  wake,  a couple  of  observations  are  in  order 
here  with  respect  to  the  applicability  of  this  recirculation  approach. 
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The  majority  of  true  gaseous  phase  stability  limit  data  may  he  correlated  with  a 
dimensional  grouping  of  the  form: 

— (f'a).  <57» 

P.NT.** 

An  example  from  Reference  12  Ls  shown  in  Figure  27. 

For  ga.seous  stability  analysis,  this  curve  follows  a fuel-air  function  of  form  similar  to  the 
kinetic  rate  function  vs  fuel-air  ratio.  Hence,  the  limit  represents  a balance  between  the 
recirculation  rate  and  the  kinetic  mass  conversion  rate.  This  then  implies  that  the  gnwiping 
Vpim/pAT.’  * represents  a grraiping  for  the  recirculation  rate. 

From  kinetics,  the  reaction  curve  follows  a shape  related  to  A/V.,p*  for  second-order 
reactions. 

For  fixed  geometry  and  flow  conditions: 


A = w. 


V„  = C.lL'DtfB'DlN* 


Combining  these  we  have: 


V„p*  " RT. 


r N p.* 


V„p*  ~ rR  pJS 

If  we  introduce  the  temperature  dependence  of  r'  as  r'  aT.".  we  obtain; 


constant  x 


IVNT.‘* 


This  is  extremely  close  to  the  data  conation  and  substantiates  the  functional  form  of  the 
recirculation  anaivsis. 


(5)  R0elrcul»tlon  Zonm  Raactfon  Alnatfet 

The  recirculation  /one  wake  of  the  bluff  body  stabilirer  is  analyzed  by  assuming  that  it 
behaves  similarly  to  a well-stirred  reactor  writh  the  volume  and  mass  entry  rate  known  from  the 
results  of  the  RECIRC  analyses. 

Analyses  of  well-stirred  reactors  have  been  presented  by  numerous  investigators  and  the 
analogy  to  actual  combustors  and  bluff  body  wakes  studied  (References  I.*!  through  211.  The  basic 
thesis  of  these  studies  is  that  the  performance  of  the  reactor  may  be  evaluated  from  a balance 
between  the  mass  entry  rate  and  the  kinetic  conversion  rate.  For  the  purpose  of  these  studies,  we 
assume  that  the  entire  wake  volume  is  available  for  reaction  and  that  the  mixing  is  very  rapid. 
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Figure  27.  Baffle  Stability  Correction 
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The?«p  assumptions  rrsult  in  an  optimistic  evaluation  of  the  kinetic  limits  which  we  will  correct 
later. 


The  following  analysis  follows  the  development  of  Kretschmer  and  Cdgers  (Reference  21 ) for 
lean  wake  fuel-air  ratios.  The  general  form  of  the  reaction  of  hydrocarbon  fuel.s  pn>ceeds  in 
essentially  two  major  steps.  The  first  is  the  pyrolysis  and  partial  oxidation  of  the  virgin  fuel  to 
form  short-lived  intermediate  species.  At  the  end  of  this  stage,  the  maximum  concentration  of 
carbon  monoxide  is  present.  The  second  stage  is  the  oxidation  of  carbon  monoxide  to  form  carbon 
dioxide.  This  latter  ste^  i,-.  much  slower  and  serves  as  the  rate-controlling  process.  The  slower  CO 
oxidation  results  in  the  very  close  similarity  of  reaction  rates  for  a wide  range  of  saturated 
hydrocarbons. 

The  general  form  of  the  conversion  equation  may  be  written  as: 


dm,,  k V B 

for  gaseous  flow  into  a well-stirred  reactor,  this  becomes: 
A k(m-H)  . ...  c" 


V.,p" 


k(m-Kl)  , 

R,"y«  T"-** 


(64) 


(65) 


For  the  assumed  single-step  reaction  process  postulated  here,  the  reaction  mass  balance  is  (for 
propane  fuel): 


0CtH4  -f  50|  -f  5mN|  — • 3yt  COt  + 4y«HjO  -(■ 
(^-y»)rjf,  -t-  5(l-y«)  O,  -f  5mN,. 


(66) 


Also,  a linear  efficiency  vs  temperature  function  is  assumed: 


T = T,  + «AT.  (67) 

Fnnn  these  equations,  the  stirred  reactor  loading  capability  may  be  written  as: 

_A_  . k(m^l)[5(l-y0j« 

V..P"  R,"y.  (5(m-H) -t  « + y«jMT,+fAT]"-**' 


Based  on  comparison  of  predicted  results  with  available  stirred  reactor  data.  Reference  21 
recommends  the  following  values  for  this  reaction: 


n:  for0<l,  n»'20 
for  ^ >1,  n=2/0 
a:  o=n/2 

C:  C ••  EyR,  See  Figure  28. 

1 

This  yields: 


np- 

for  lean  mixtures. 


1,29  X 10‘»  (m  + 1)  f5(l-yc)|*  (<»-yt)*  e 
(0.08206)'*  y«(6(m-H)  -f  0 -t  ytp  (T.-fi AT| 


(69) 


sHwm 
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Figure  28.  Variation  in  Activation  Energy  with  Inlet  Temperature  and  Equivalence  Ratio 

Wliilp  tlip  results  presented  above  are  from  one  reference,  siniilar  results  have  been  obtained 
by  others  tor  the  same  problem. 

lion^twell  and  Weiss  (Heferenee  M)  present  the  following  results  for  lean  and  rieh  operation: 


[.can: 


V„p'>  ' Tr”  <ti<|4.7f)  + - <)]■• 

Rieh: 

_A_  . 4.,o  k ^ 

V„p"  ' ' 'IV  ’ ( I 4.76  - 1 + 0.080  (1  + 1G() 


(70) 


(71) 


where 


ki  - 1.67  X 10'°  litro”  "/‘’K”  ‘ !;m-mole“*  sec 
kj  = l.n  X 10”  litre‘'  V°K'’‘ f;m-mole‘”,sec 
M ~ 42,0(K)  );m-cal/(t  n-mole. 

'I  he  solution  iiroeedure  utili?;ed  in  the  WAKK  subrotitine  utilizes  the  results  of  RKCIRC  as 
a definition  of  the  entry  rate,  A,  and  the  zone  Volume,  V,„  the  eriuivalenee  ratio  is  from  Equation 
20  and  the  re.sulls  of  111,  112, 116  and  RECIRC.  The  ideal  temi)crature  rise,  AT,  is  read  from  curves 
of  AT  versus  wake  fuel-air  ratio,  inlet  temperature,  and  .static  pressure  for  the  specified  fuel  type. 
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An  iterative  solution  procedure  solves  for  the  efficiency, «,  where  the  known  value  of  AA^„p*  agrees 
with  the  predicted  value  as  a function  of  f.  The  complexity  of  Equation  70  or  71  requires  this  sort 
(rf  reverse  solution. 


The  behavior  of  the  solution  is  shown  in  Figure  29  for  inlet  temperatures  of  300  and  400®K 
and  equivalence  ratios  of  0.8  and  1.0.  The  wake  efficiency  decreases  at  an  increasing  rate  until  the 
decay  slope  becomes  infinite.  At  this  point,  the  wake  reaction  pitxress  is  said  to  have  blownout. 
If  we  plot  the  locus  of  the  blowout  points  as  A/VoP*  versus  equivalence  ratio,  we  obtain  a clas.sical 
blowout  curve.  This  curve,  except  for  magnitude,  repr3sents  the  classical  results  for  gaseous  fuel 
data  from  wake-stabilized  flames.  This  is  shown  in  Figure  30. 


This  analysis  tends  to  overestimate  the  limits  of  blowout  velocity  when  compared  to 
available  data.  To  reconcile  this,  comparisons  were  made  between  predicted  limits  and  actual 
limits  for  available  data  (References  22  through  28).  The  correlating  ratio  between  predicted  and 
actual  was  3.55,  i.e.: 


( vt*  ) 


mux.  prvd 


i^) 


— = 3.56. 


max 


(72) 


In  the  analysis,  the  calculated  recirculation  rate  per  unit  volume  is  multiplied  by  this  factor 
before  solving  for  t.  In  this  manner,  the  predicted  limits  and  available  data  are  numerically  equal, 
and  the  efficiency  response  towards  blowout  follows  the  theoretical  curve. 

The  reason  for  this  error  is  most  likely  either  imperfect  mixing  or  lack  of  full  utilization  of 
the  wake  volume  for  reaction.  The  approach  used  above  will  account  for  either  one  of  these. 


Figure  29.  Stirred  Reactor  Kinetic  Efficiency  vs  Loading  Rate 
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(6j  Turbulent  Flame  Spreading  Rate 

The  turbulent  flame  propagation  into  the  unreacted  free-stream  is  initiated  in  the  shear 
layers  of  the  wake.  The  model  used  (subroutine  ST)  related  the  local  turbulent  flame  speed  to  the 
local  aerothermodynamic  conditions  and  performs  a finite  difference  integration  of  the  flame 
front  penetration  starting  in  the  wake  and  proceeding  to  the  exhaust  nozzle. 

For  the  purposes  of  current  analysis,  the  following  assumptions  were  made; 

1.  Uniform  airflow  profiles 

2.  Uniform  fuel-air  ratio 

3.  Incompressible  acceleration  of  free  air  velocity  by  the  flameholder  blockage 
with  no  induced  profile 

4.  Known  wake  size  and  reaction  efliciency 

5.  Two-dimensional  ducted  flame. 


The  schematic  of  the  situation  which  is  analyzed  is  shown  in  Figure  31. 


The  approach  flow,  at  known  levels  of  pressure,  temperature,  velocity  and  fuel-air  ratio,  is 
accelerated  by  the  blockage  of  the  flameholder  to  velocity  U,  where; 


U 


V. 

(l-D 


where: 


(73) 


U ~ Velocity  at  flameholder  tip 
V,  -»  Approach  velocity 
r ~ Blockage  ratio. 


At  this  point.  Station  1 of  Figure  31,  an  induced  turbulence  level  is  calculated  from: 

)■}!]” 


(74) 


This  equation  (Reference  12)  relates  the  turbulence  intensity,  to,  to  the  blockage  ratio  and  the 
flameholder  zero  blockage  drag  coefficient,  C^. 

At  this  location,  the  turbulent  flame  velocity  calculations  are  initiated.  The  equation  used 
for  the  local  flame  speed  is  the  Karlovitz  equation  (Reference  22): 

St  = S,  (2u'S  )*•  (75) 

where: 

St  ~ Turbulent  flame  speed,  ft/sec 
Si  ~ Laminar  flame  sp^,  ft/sec 
u'  ~ RMS  turbulence  velocity,  ft/sec. 


\ 
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Air  Velocity  Accelerated 
By  Blockage  - Uniform 
Profile  Assumed 


Figure  31.  Schematic  of  Flame  Spreading  Analysis 


The  value  of  u'  is: 


u'  = «„  U. 

An  additional  term  is  required  to  relate  the  resultant  flame  speed  to  the  level  of 
recirculation  zone  reaction  efficiency.  As  the  wake  efficiency  and  temperature  decrease,  local 
areas  in  the  flame  sheet  appear  where  ignition  does  not  occur.  Since  the  model  assumes  a 
continuous  flame  sheet,  we  model  these  local  areas  of  no  ignition  by  reducing  the  overall  flame 
speed.  A reduced  f.ame  speed  results  in  lower  overall  combustion  efficiency  wbicb  is  similar  to  tbe 
effect  of  locally  zero  ignition.  The  model  uses  the  following  equation: 

Sf  = St  • i/w.  (77) 

The  initial  value  for  the  augmentor  efficiency  is  the  wake  reaction  level  on  a mass  weighted 
basis.  Expressed  as  an  equation  this  is: 

>?.•„  I?*  ' ( <8) 

m, 

wiiere: 

n,.  ~ Initial  efficiencv 

t) 

ri„  ~ Wake  efficiency 
m,  ~ Wake  mass  flowrate 
m,  ~ Total  duct  flowrate. 

The  type  of  flame  utilized  in  this  model  is  a zero  thickne.ss  flame  which  separates  a region 
of  unreacted  propellants  from  a region  of  completely  reacted  products.  From  this  set  up  the 
average  local  augmentor  efficiency  is  simply  the  ratio  of  the  transverse  flame  penetration,  Ay,  to 
the  duct  width,  W,  (see  Figure  .31). 

To  be  consistent,  the  transverse  location  of  the  flame  front  at  the  initial  calculation  station 
is  taken  to  be: 

Ay.,  = ■ W (79) 

This  value  is  assigned  to  the  first  axial  .station,  which  is  assumed  to  occur  halfway  down  the 
length  of  the  recirculation  zone.  From  visual  observations  of  wake-.stabilized  flames,  this  is  the 
approximate  location  of  transverse  flame  initiation. 

From  this  location  dowi  stream  to  t.ie  exhau.st  nozzle,  the  flame  front  transverse  location  is 
calculated  by  a finite  difference  integration  of  the  local  flame  speed.  Several  axial  profiles  are 
introduced  as  the  integration  proceeds.  These  are: 

• The  turbulence  intensity  is  decayed  from  the  value  generated  at  the  aft 
flameholder  lip  (Equation  74)  at  a rate  inversely  proportional  to  the  square 
root  of  axial  distance  over  an  effective  jet  length.  The  final  value  is  set  at  the 
initial  turbulence  level.  The  effective  jet  length  is  set  at  10  IVD  where  the  D 
is  the  open  area  distance  between  adjacent  flameholders. 

® The  velocity  of  the  nnreacted  fuel-air  mixture  is  retained  at  the  level 
generated  at  the  flameholder  lip.  Measured  profiles  from  several  ducted 
flame  test  rigs  support  this  a.ssumption. 
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• A tenn  is  introduced  which  relates  the  local  flame  speed  to  the  local  average 
duct  combustion  efHciency,  peaking  at  50%.  This  treats  the  counteracting 
influences  of  reduced  heat  loss  as  efficiency  increases  and  the  reduced  free 
oxygen  concentration.  Local  rates  which  follow  roughly  a sine  wave  function 
have  been  reported  from  duct  data. 

An  additional  term  is  added  to  account  for  the  reduction  in  flame  speed  of  a fuel  spray 
compared  to  a premixed  flame.  This  term  relates  the  ratio  of  eflective  flame  speed  to  premixeu 
laminar  flame  speed  following  the  results  of  Reference  30.  It  accounts  for  the  complicated 
interactions  during  flame  spreading  in  an  evaporating  spray  in  a simplified  manner.  The  effect  of 
the  liquid  droplet  diameter  is  shown  in  Figure  32.  The  droplet  diameter  utilized  in  the  analysis 
will  be  the  mean  diameter  as  it  exists  at  the  flameholder  trailing  edge. 


Figure  32.  Flame  Speed  for  Monodisperse  Tetralin  Spray 


Analysis  of  the  terms  utilized  for  evaluation  of  the  laminar  flame  speed  term,  S'„  has 
res^.ved  in  the  following  (Reference  23): 

(■&)“  (^)'  '«’• 


where: 

$ = laminar  flame  speed  at  1 atm  and  540° 

0 = equivalence  ratio 
T,  = air  temperature,  °F 
Xo,  = oxygen  mole  fraction. 

The  influence  of  pressure  ih  indeterminate  at  this  time  and  has  been  incorporated  as  s/p  for 
subatmospheric  data  and  no  influence  for  pressures  above  1 atmosphere. 


( 
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Several  results  have  been  generated  from  this  model  in  its  current  stage  of  evaluation.  The 
first  is  shown  in  Figure  33  for  a 60-in.  duct,  2-in.  flameholder,  and  2.5'';  blockage  at  2.50  ft/sec 
entrance  velocity.  The  curve  shows  a significant  rate  of  increase  in  efficiency  between  0.0.50  and 
0.060  fuel-air  ratio  and  a steep  decline  past  0.090  f/a. 

Figure  .34  shows  the  local  values  of  efficiency  versus  axial  length  for  a fixed  duct  width  and 
various  flameholder  widths.  For  this  case,  the  blockage  is  proportional  to  the  width.  The 
interacting  effects  of  increased  lip  velocity  and  increa.sed  induced  turbulence  are  such  that  the 
net  effect  is  reduced  overall  efficiency  as  the  blockage  is  increased.  Since  in  this  case  complete 
wake  reaction  efficiency  is  assumed,  somj  of  this  effect  will  be  reduced  by  the  decreased  wake 
efficiency  as  the  flameholder  size  is  reduced.  This  blockage  effect  is  in  agreement  with  results 
from  the  experimental  portion  of  this  program. 

The  effect  of  free-stream  turbulence  is  shown  in  Figure  35.  The  anticipated  effect  of 
increased  flame  speed  generated  by  the  higher  turbulence  was  verified. 

The  effect  of  reduced  wake  reaction  level  on  the  overall  efficiency  is  shown  in  Figure  .36.  The 
relationship  is  such  that  near  the  peak  flame  speed  le\els,  a unit  reduction  in  wake  efficiency 
results  in  less  than  a unit  reduction  in  overall  efficiencj’.  This  implies  a reduced  sensitivity  of  the 
overall  efficiency  to  oscillations  of  the  wake  efficiency  relative  to  what  was  previously  assumed. 

A companion  prosram  to  this.  Contract  F.3.361.5-76-C-2024.  utilized  a correlation  for 
augmentor  efficiency  called  the  “Beta  Correlation.”  Since  this  correlation  represents  a fairly  large 
body  of  teat  data,  comparison  of  the  results  of  this  subroutine  with  the  Beta  Correlation  is  a quick 
way  of  comparing  the  mc.del  to  test  results.  Two  comparisons  are  shown  in  Figure  .37.  As  shown, 
good  agreement  is  achieved. 


Fuel-Air  Ratio 

ro  i.wBs 


Figure  33.  Average  Efficiency  vs  Fuel  Air  Ratio 
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Figure  36.  Effect  of  Wake  Reaction  Level  on  Duct 
Combuefion  Level 
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2.  EXPERIMENTAL  INVESTIGATION 


a.  Test  Apparatus 

Tht  study  experiments  were  conducted  in  a (boiler  plate)  comb'i»rii>n  sv>.tem  in  whii  h fh** 
pressure,  temperature,  jtas  flow  and  flameholder  instability  in  an  aiurmentor  were  stmul.itefl  The 
circular  cross-sectional  rig  was  designed  and  built  under  an  indefiendent  f*ratt  A-  \Vhitne\ 
Aircraft  Group  research  program.  The  simulator  was  designed  .<o  that  vahisis  I'onfig’i'ari'in 
chariges  including  flameholders.  spraybars.  and  distances  between  reflective  points  couhl  in* 
easily  made  f«)r  diagmwtic  evaluation 

A drawing  of  the  test  rig  Ls  shown  in  Figure  The  rig  is  made  up  of  an  inlet  refle»  t ive  onfi,  e 
plate,  inlet  ca.se.  fuel  injection  case,  combustion  case,  transitii'n  flange,  and  a variable  area 
exhaast  nozzle.  Photographs  of  the  test  hardware  are  shown  in  Figures  thnnigh  47. 

The  orifice  plate  (Figure  ."fHl  pnnides  the  system  upstearn  reflective  point  The  overall 
system  length  can  be  varied  from  Ib'J  4 to  1 19.4  inches  by  the  location  chosen  for  the  orifice  plate 
It  can  be  placed  at  the  entrance  or  exit  of  the  inlet  case.  The  inlet  case  i Figure  4<)i  is  :{9  9 ini  hts. 
long.  When  inlet  temperatures  above  ftiO’F  are  desired,  the  inlet  case  is  repla«-ed  with  a heater 
burner  that  consists  of  a burner  can  and  fuel  nozzle  cluster  mounted  in  a fioiler  plate  hiKising  .All 
high  temperature  (W**  to  l.'KiO'Fi  testing  was  done  with  a system  length  of  1 19  4 inches. 

The  uncooled  fuel  injection  section  (Figure  41f  is  99  4 inches  li>ng  arid  has  two  tntecfion 
planes  which  can  he  fitted  with  either  three  liquid  fuel  *>r  three  gaseous  fuel  spravbars  Tvpi«  al 
liquid  and  gaseous  spraybars  are  shown  in  Figure  42.  The  primarv  and  sei-ondarv  injei  tion  planev 
are  knated  H and  40  inches,  respectively,  upstream  of  the  pnmarv  flameholciee  Im  arion  The 
zone  t spraybar  is  located  on  the  rig  centerline.  The  zone  2 and  zone  9 spravbars  are  located  2.9 
inches  outboard  of  the  rig  centerline.  The  location  and  the  spray  direction  of  the  liquid  fuel  zones 
are  .shown  m Figure  4;1.  The  injection  direction  for  each  zone  was  the  same  for  the  ga.-  sprav  bars 

The  comhu-stion  section  (Figure  441  is  a water  cooled  duct  12  inches  in  diai.seter.  T*i  0 inches 
long  and  has  a primary  and  secondary  flameholder  kxation.  The  primary  flat.ieholder  loi  ation  is 
79.5  inches  from  the  nozzle  exit  plane  and  th'  secondary  location  is  47. .5  inches  from  the  no/zle 
exit  plane.  Either  position  can  he  fitted  writh  three  “V"  gutter  flameholders  of  .9.5  and  52'. 
blockage,  1:1  and  2:1  length  to  width,  and  draft  angles  of  22.5  and  45'  An  oxy-acetylene  pilot 
burner  is  mounted  in  the  recirculation  zone  of  the  center  flameholder.  This  burner  ran  pro\  ide  a 
rontinuoas  ignition  source  and  was  used  to  supply  the  variations  in  heat  addition  The  various 
flameholders  configuration.^  are  shown  in  Figure  45.  The  centerline  of  the  three  flameholder  zones 
are  aligned  with  the  centerline  of  the  respective  fuel  spraybar  zones. 

The  water  cooled  transition  flange  (Figure  46*  ■ hanges  the  rig  flowpath  fn»m  circular  to 
rectangular  to  match  the  variable  area  nozzle.  The  exl-a  ist  nozzle  assembly  is  designed  to  permit 
cnntinunas  variation  of  exit  area  over  a wide  range  «h  choked  operating  conditions.  The  nozzle 
assembly  consists  of  (1)  two  remotely  actuated  tylmdrical  water-cooled  rods.  4'..  inches  iu 
diameter  on  one  end  and  I ’l  inches  in  diameter  on  the  other.  (2>  two  semicylindricat  side  wall 
plugs,  and  (3)  two  flat  plate  sidewalls.  Moving  the  rod«  in  and  out  of  the  duct,  with  and  withcxit 
the  side  wall  plug>,  results  in  a geometric  area  change  ranging  fn>m  16. 1 to  1 12.9  in.*  Figure  47 
shows  the  possible  geometric  area  extremes.  The  rods  have  been  provided  with  total  pres.sure 
ports  on  the  upstream  side  and  static  pressure  taps  on  the  sides  of  the  rods,  so  that  at  any  rod 
position  there  is  a minimum  of  five  total  pressure  and  four  static  pressure  pickups  per  rod  in  the 
duct.  As  shown  in  Figure  48.  each  rod  is  driven  independently  by  a linear  actuator  powered  by  a 
24  volt  dc  motor. 
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•e  39.  Rig  Hardware 


Direction  and  Pintel  Location  with  Liquid  Fuel  Injectors 


Figure  44.  Combustion  Cnr; 


Figure  45.  Flameholders 


Figure  47.  Geometric  Nome  Area  Extremes 


b.  Nozzle  Actuator  Assembly 


b.  Instrumentation 


Thr  ifunrumentation.  iUled  in  Tahie  2 u-^  in  in<mitr>r  the  foUoskini;  tr<>t  parameter* 

• Airflow 

• Fuel  flow 

• Ri):  inlet  total  pre>4ure  and  temperature 

• CombuMfir  inlet  total  preiwure  and  temperature 

• romb-.tst.-s’  exit  total  preswure 

• Comhuntor  pressure  iKcillations 

• F’-ei  temperature,  flowrate,  and  pressure 

• Velos  itv  changes  in  duct 

• Temfierature  change*  in  comhuslor 

• Flameholder  skin  temperatures 

• l*il«»f  burner  ox\gen  and  acetviene  fl«*w  rates 

3 Fuel  ai*  ratio  upstream  ot  the  flameholder  and  in  the  flameKilder  wake 

The  data  aere  re<-orded  bv  an  automatic  data  rect»fding  *v*tem  and  redu«ed  through  the 
IBM  :I70.  Ifvs  computer  *\stem 

Special  test  item*  unique  to  this  program  are  discussed  in  the  following  paragraph* 

The  combustK  was  instrumented  with  fi%e  high  resp else  Ki*tler  M<del  fiUKA  presauie 
transducers  (Fig'jre  •I9i  to  determine  hrgh  the  frequency  and  type  of  wave  «iccumng  dunng 
combustion  instahi.'ity  The  axial  and  angular  locations  for  each  of  the  Kistlers  are  sh>"«m  in 
Figure  ^ The  phase  relationship  and  relative  amplitude  of  the  iwcillatKins.  sensed  hy  the 
transd'jcers  brated  at  intervals  down  the  hg.  provided  the  necessary  information  to  identify  the 
wave  pattern  and  amplitude  gains  dunng  rumble  operaton 

An  oxv  aietviens-  pilot  tiech  wa*  mounted  in  the  rerimilatioo  rone  i4  the  renter 
flameholder  This  uevh  provided  continuous  ignitioo  and  •imulafed  piloting  hv  wake  heat 
addition 

The  center  flamelV'Ider  pwilion.  rone  1.  was  instrumented  with  tsro  skin  thrrm'>ro«iples, 
Thev  were  attached  io  the  upstream  side  of  the  "V  gutter  on  the  I M • 1 .2.  1 .2  ■ 1.2.  and  1 .2  • 
0.6  flameholder  (onriguraiams  A typjca!  thermocouple  location  is  shosm  in  Figure  .V», 

Two  four-inch  diameter  vveor  gla*e  viewing  poets  swe  located  !*>'  apart  at  the  trailing 
edge  of  the  pnmarv  flameholder  so  that  high  speed  coi<e  m<aH«)  pclure*  could  he  made  dunng 
sieadv -slate  and  rumble  operateets  to  pmside  comparative  data  on  flow  dvnamic-  and 
roaihustain  Two  air.c<a>led  Keuratmn  pnihes  (Figure  Wi  are  u«ed  to  determine  ga*  temperature 
increase  •»  decrease  The  temperature  changes  are  correlated  with  vebeitv  changes  detected  hy 
a strain  gage  attached  ft.  a target  probe  upstream  of  the  flameholders  The  oereiatMin  between 
Inca'  temperature  and  velocitv  changes  provide  additional  data  lei  the  rumble  mechanism. 

The  fuel  vapeirataei  probe  »%  located  |usl  upstream  <if  the  flameholder  and  ran  he  tnive»sed 
actiMS  the  flameholder  • Figure  .VI  The  purpose  «if  this  probe  is  to  sample  the  fuel -air  mixture  just 
upstream  of  the  flameholder  and  provide  a mewns  of  quantifying  the  aresaint  of  liquid  and 
vapaired  furl  present 
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Sides  of  the  Water-Cooled  R(M>m  MonitorioK.  Automalii- 
Nozzle  Rods.  (Six  in  each  Data  Recording. 


Figure  50.  Flameholder  Skin  Thermocouple  Locations 


c.  Test  Program 

Concurrent  with  Task  I analytical  studies,  an  experimental  program  was  conducted  to 
verify  the  analytical  model  and  to  determine  values  for  portions  of  the  flame  stabilization  model. 
Studies  were  conducted  over  a range  of  fuel-air  ratios  to  evaluate  the  influence  of  gutter  geometry 
and  blockage,  pressure,  temperature,  velocity,  and  fuel  vaporization  level.  Additional  tests  were 
made  to  evaluate  the  influence  of  heat  addition  to  the  wake  and  the  influence  of  inclined 
flameholder  gutters. 

Planned  test  series  (Figure  51)  were  designed  to  cover  a range  of  inlet  conditions  typical  of 
augmentors,  especially  the  conditions  which  simulate  the  effect  of  the  cold,  fan  duct  airstreani  in 
turbofan  engines.  Particular  emphasis  was  placed  on  the  cold  air  influence,  since  turbofan 
augmentor  experience  at  Pratt  & Whitney  Aircraft  has  shown  that  rumble  may  be  reduced  or 
eliminated  by  alterations  to  the  cold-stream  flameholder  alone. 

The  test  matrix  included  varying  the  fuel-air  ratio,  blockage  ratio,  flameholder  geometry, 
inlet  temperature,  pressure,  fuel  distribution  and  vaporization.  For  each  sequence  of  test 
variables,  a sampling  probe  was  installed  to  measure  the  local  total  and  vapor  phase  fuel-air 
ratios  upstream  of  the  flameholder.  For  several  tests,  two  sampling  tubes  were  mounted  in  the 
downstream  wall  of  the  center  flameholder  for  measurements  inside  the  near  wake  region.  A 
Beckman  Model  402  hydrocarbon  analyzer  was  used  to  measure  on  line  the  data  from  the 
upstream  probe  and  the  wake  tubes. 

After  completion  of  the  first  test  seres,  the  test  program  was  modified.  I’he  revised  test 
program  shown  in  Figure  52  was  a composite  program  combining  the  test  requirements  of  both 
the  Flameholder  Combustion  Instability  Study  and  the  companion  system  model,  Contract 
F33615-76-C-2024,  into  a single  test  program  to  cost  effectively  gather  the  most  data.  The  test 
matrix  was  not  significantly  changed  from  the  original  program,  i.e.,  most  of  the  original  items 
to  be  investigated  were  still  included.  The  only  item  deleted  was  the  effect  of  fuel  injection  system 
stiffness  because  (1)  a concensus  of  the  investigators  believed  its  impact  on  the  program  was 
minimal,  and  (2)  more  test  time  could  be  allocated  to  the  other  suspected  mechanisms.  Most  of 
the  test  program  changes  were  reflected  in  the  test  operating  conditions.  In  the  original  test 
matrix,  bands  of  pressure,  temperatures,  duct  Mach  numbers,  and  fuel-air  ratios  were  given.  The 
Series  I test  results,  however,  specifically  defined  the  test  conditions  of  interest  within  the 
capabilities  of  the  system.  The  resultant  test  conditions  are  shown  in  Table  3. 

d.  Experimental  Program  Problems 
(1)  Series  I Totting 

The  primary  purpose  of  the  Series  I testing  was  to  determine  stability  characteristics  of  the 
test  rig,  which  would  serve  as  a baseline  for  the  Series  II  test  program.  A secondary  purpo.se  of  the 
Series  I testing  was  to  evaluate  the  test  rig  hardware  and  instrumentation  which,  for  the  most 
part,  had  not  undergone  previous  testing. 
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Figure  52.  Flameholder  Model  Rig  Tests 

i 

TABLE  3 

REVISED  TEST  CONDITIONS 

Inlet  Iitlet  Equivalence  Duct 

Test  Pressure  Temperature  Ratio  Mach 

Number  (psia) (^F) 07] Number 

A 10  *200  0.5,  1.0,  1.6  0.088 

B 10  200  0.6,  1.0,  1.6  0.185 

C 10  *400  0.5,  1.0,  1.6  0.088 

D 10  400  0.6,  1.0,  1.5  0.155 

E 16  *400  0.5,  1.0,  1.6  0.088 

F 15  400  0.6,  1.0,  1.6  0.165 

G 16  400  0.6,  1.0,  1.5  0.238 

* Riin  2 or  3 torch  flowrate  variationa  at  the  completion  of  the  iao-kinetic  and 
wake  (/s  data  acquisition. 
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The  overall  test  ri;:  eoridition  after  Series  I testing,  was  very  gwid.  No  major  damage  was 
incurred.  The  target  probes,  ionisation  probes,  viewing  window,  and  a total  pressure  probe, 
however,  were  damaged  during  the  course  of  the  testing.  The  damaged  probes  were  repaired  with 
slight  modifications  incorporated  to  improve  durability  before  additional  testing.  The  set  of  vvcor 
glass  viewing  windows  was  destroyed  early  in  the  first  series  of  tests.  It  was  not  replaced  during 
the  remainder  of  the  Series  I program  because  more  valuable  visual  data  could  be  obtained  during 
later  testing  after  the  stability  characteristics  of  the  test  rig  were  known.  Another  problem  area 
uncovered  during  the  Series  I testing  was  the  upstream  orifice.  The  orifice,  which  had  been 
installed  to  provide  a known  reflective  location,  was  sending  a jet  of  airflow  throughout  the  burner 
.system  at  high  airflow  rates  ( 10  to  l.T  Hbm/sec).  This  jet  was  apparently  not  attached  to  any  wall 
surface  for  the  entire  length  o*  the  burner  -system.  To  remedy  the  problem,  prior  to  Series  II 
testing,  the  orifice  plate  was  reworked  by  plugging  the  large  single  hole  and  rematching  to  include 
sixteen  ( Hi)  smaller  holes  to  maintain  the  same  open  area.  Since  combustion  efficiency  could  not 
be  determined  from  the  data  with  an  unchoked  nozzle,  the  sidewall  plugs,  di.scussed  earlier 
(Figure  47),  were  installed  in  the  nozzle  before  Series  II  testing  to  ensure  that  choked  flow  could 
he  maintained  and  combustion  efficiency  determined  over  the  full  range  of  operating  conditions. 
The  two-phase  sampling  probe  was  not  operational  during  the  first  series  of  testing. 

{2)  Series  11  Testing 

Some  minor  problems  were  also  encountered  during  the  Series  II  testing.  The  test  facility 
ejector  .system  was  not  operating  at  specified  efficiency  and  would  have  required  major  repairs  to 
correct  the  deficiencies.  Since  combustion  efficiency  measurements  were  desired,  choked  flow  at 
the  exhaust  nozzle  was  required.  To  maintain  choked  flow,  the  lowest  rig  pressure  was  limited  to 
approximately  14  psia  and  the  rig  duct  Mach  number  to  0.1 2 compared  to  a planned  rig  pressure 
of  10  psia  and  duct  Mach  number  of  O.l.'i.S.  It  was  determined  that  these  revised  conditions  would 
provide  the  required  input  for  the  .aodel  and  the  te.st  program  was  modified  accordingly. 

Upon  completion  of  the  test  program,  it  was  determined  that  the  target  probe,  which  had 
been  modified  at  the  onclusion  of  the  Scries  I te.sting,  was  often  vibrating  at  frequencies  between 
200  and  6tK)  Hz  when  the  airflow  dynamic  pressure  measurements  were  55  Hz.  It  was  also 
determined  that  the  ionization  probes  were  occasionally  grounded  due  to  dampness  that  resulted 
in  a 60-cycle,  high-gain  signal  Because  the  te.st  program  had  been  completed  when  these  two 
discrepancies  were  determined,  and  the  elimination  of  these  two  pieces  of  instrumentation  would 
have  only  a minor  impact  on  the  model  input,  the  data  was  disregarded.  It  was  felt  that 
disregarding  all  of  the  data  would  be  more  prudent  than  drawing  conclusions  from  data  that  may 
or  may  not  be  correct. 

The  final  problem  area  of  Series  II  testing  involved  the  zone  3 spraybar.  After  completion  of 
the  testing  with  turbulence  screens,  it  was  determined  that  the  spraybar  had  been  partially 
plugged  with  teflon  tape  that  had  been  used  to  seal  the  fuel  .system  fittings.  The  tape  was 
removed  and  the  spraybar  recalibrated  prior  to  resuming  the  te.st  program,  but  the  data  from  the 
previous  test  sequence  involving  zone  3 fuel  flow  was  disregarded.  This  problem  also  had  a minor 
impact  on  the  model  input. 


3.  EXPERI.MEMTAL  PROGHAM  RESULTS 
a.  Ovornll  Commonts  and  Roaulta 

The  experimental  program  was  structured  to  evaluate  the  impact  on  the  combustion  process 
of  the  following  major  variables: 

o Fuel-air  ratio 
a Spray  vaporization 
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• Fuel  distribution 

• Heat  addition  to  the  wake 

• Combustion  length 

• Inlet  conditions,  pressure,  temperature,  velocity 

• Flameholder  geometry  and  blockage 

• Inclined  flameholders. 

The  test  matrices  presented  in  Section  11-2  included  a fuel-air  ratio  excursion  at  several  pressures, 
temperatures  and  flow  velocities  for  each  majo-  configuration  change.  Several  of  these  were  run 
with  and  then  without  an  oxyacetylene  torch  m mted  to  provide  heat  addition  to  the  wake  of  the 
central  flameholder. 

The  major  results  of  the  test  program  may  be  summarized  as  follows: 

• It  was  not  possible  to  maintain  a stable  combustion  process  in  the  system  at 
approach  flow  Mach  numbers  above  approximately  0.12  at  inlet  tem- 
peratures of  200  to  400°  F. 

• The  teat  data  at  the.se  lower  temperatures  exhibited  a large  scatter  and  a lack 
of  definitive  repeatability  with  regard  to  efficiency.  As  such,  the  recorded 
efficiencies  should  be  interpreted  as  accurate  only  to  about  ± lO'^f  absolute. 

• A definite  correlation  was  observed  between  the  recorded  metal  temperature 
of  the  flameholder  and  the  overall  combustion  efficiency.  The  correlation 
varied  from  one  system  configuration  to  another  but  was  fairly  good  for  each 
system. 

• The  onset  of  rumble  and  the  decline  in  combustion  efficiency  occurred  at  the 
same  time,  reinforcing  the  combustion  efficiency-driven  oscillation  concept. 

Before  going  further  into  the  data  analysis  a test  case  will  be  examined  to  substantiate  the 
applicability  of  the  two-phase  fuel  model  concept. 

The  data  here  are  from  the  first  run  of  Series  II  after  modii't  utions  to  the  variable  area 
exhaust  nozzle.  This  run,  7.01,  utilized  a 52%  blockage  flameholder  w ith  1.8-in.  wide  gutters.  The 
sepcra^'oii  between  flameholders  and  spraybars  wes  8 in.  The  me:‘sured  overall  combustion 
eff  -iency  j.d  r-t  /.ble  amplitude  are  shown  verses  fuel-air  ratio  in  Fij  ;res  53  and  54.  The  results 
indicate  s'gnificanv  ..aviations  from  premixed  behavior  in  that  peak  eff.  'iencies  occur  at  less  than 
p' oichiometric  fwl-air  ratios. 


The  two-phase  fuel  wake  composition  equation: 

I"- /?,  + (l-(3.)-^  (81) 

allows  for  the  existence  of  a wake  vapor  fuel-air  ratio  which  is  on  the  rich  side  of  stoichiometric 
( 0 > 1.0)  even  when  the  injected  fuel-air  ratio  is  lean  (0  < 1.0).  As  the  injected  fuel-air  ratio  was 
increased,  the  wake  also  increased  and,  being  rich,  the  wake  reaction  efilciency  decreased.  Thus, 
the  overall  efficiency  would  exhibit  a peak  value  while  the  injected  fuel-air  ratio  was  lean.  For  the 
conditions  of  Run  7.01  (i.e.,  400°F  inlet  temperature,  near  ambient  pressure  and  150  fps  velocity) 
the  calculated  value  of  di  is  32%.  The  droplet  averaged  collection  efficiency,  /?„  is  calculated  at 
85%.  The  wake  recirculation  rate,  K„  for  the  52%  blockage  is  23%.  The  percentage  of  the  surface 
film  which  vaporizes,  0„  is  calculated  to  be  about  40%  for  a fuel-air  ratio  of  0.050. 
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Figure  54.  Rumble  Break  Point  Occurs  at  Predicted  Point  of  Rich  Wake  Transfer 
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Factoring  these  into  equation  81  yields: 


(0.32)  -r  (0.68)  — (H2» 

or  a wake  fuel-air  ratio  of  0.066  at  the  inlet  fuel-air  ratio  of  O.O.'iO. 

As  the  overall  fuel-air  ratio  increases,  the  wake  fuel  air  ratir  o tr  '-u  ps  to  pace  it  at  32'- 
greater  value.  At  0.070  overall  fuel-air  ratio,  the  wake  is  at  0.0.52,  Due  tii  the  veiy  rapid  decrease 
in  wake  reaction  capability  as  the  fuel-air  ratio  exceeds  stoichiometric,  the  overall  .system 
efficiency  decreases  as  we  go  from  0.050  to  0.070  overall  .'uel-air  ratio. 

The  data  from  Run  7.01  indicate  peak  efliciency  at  approximately  a 0.052  f/a.  A premixed 
system  would  produce  a peak  value  at  e’  .ut  107'';  stoichiometric,  or  0.0728  f/a.  This  ratio  is  1 ..399 
and  is  very  close  to  the  predicted  v ke  enrichment  ratio  of  1.32.  These  results  confirm  the 
anticipated  results  of  the  prr  V .ions. 

■^he  influence  1 ihe  fuel  vaporization  level  on  the  efficiency  was  studied  by  two  ways.  A 
purely  gas^-  .'ovi  case  was  run  utilizing  natural  gas  (98'^5  CH,)  rather  than  the  liquid  •JP-4.  and 
a test  was  run  with  liquid  fuel  and  the  spraybar  to  ffameholder  separation  increased  from  8 to  40 
in.  This  ] .ter  shift  would  increase  the  level  of  fuel  droplet  vaporization  from  32  to  64''r  at  the 
400°F  inlet  condition.  This  shifts  the  enrichment  ratio  (equation  81)  from  1.32  to  1.17.  This  te.st 
was  Run  9.01  and  the  400°F  results  are  shown  in  Figure  55.  For  this  case  the  pronounced  peak  at 
fuel-air  ratios  below  stoichiometric  is  absent.  The  model  predicts  the  peak  efficiency  fuel-air  ratio 
to  be  0.062. 
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Figure  55.  Overall  Combustion  Efficiency  us  Fuel-Air  Ratio 
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The  longer  mixing  and  vaporization  length  test  pitidured  larger  data  scatter  and  higher 
rumble  amplitudes.  It  is  felt  that  these  results  are  due  to  the  increased  time  delay  between 
injection  and  stabilization.  The  longer  time  delay  allows  eny  gas  velocity  oscillations  to  pn>duce 
greater  oscillations  in  the  droplet  velocity  through  the  drag  terms.  These  liquid  velocity 
oscillations  generate  fuel-air  «)scillations  at  the  flameholder  which  are  ultimately  seen  as 
efficiency  oscillations,  or  rumble. 

The  largest  problem  which  exists  with  the  data  as  taken  is  the  rather  severe  fuel  distribution 
errors  which  were  pre.sent.  The  location  of  the  flameholders  and  the  flow  characteri.stics  of  the 
spraybars  resulted  in  the  center  flameholder  being  run  at  a fuel-air  ratio  several  times  greater 
than  the  outer  ones.  This  distribution  requires  that  analyses  over  the  overall  efficiency  data  be 
done  on  a multiple  streamtube  basis. 

The  flameholders  are  shown  in  Figure  o6  for  the  -Vi'v  blockage  flameholder  (1.8-in.  widel. 
The  open  area  distribution,  from  top  to  bottom,  is  l-'i.l.  11.6.  1 1.8.  I.t.1  sq  in..  As  a percentage  this 
is  28.  22,  22,  28.  For  the  three  1 .8-in.  wide  flameholder  gutters,  assuming  that  the  flow  streamline 
split  between  adjacent  gutters  on  an  equal  basis,  the  effective  blockages  are  47.6.  6-^.  47.6.  The 
airflow  distribution.s  assigned  to  the  three  gutters  becomes  39.  22  and  39  r)f  the  total  airflf)w. 

Concurrent  with  this,  the  fuel  flowrate  from  the  three  spraybarr.  was  such  that  the  centerline 
spraybar  was  producing  about  twice  the  fuel  flow  of  either  of  the  other  two  spraybars.  For  a fuel 
flowrate  distribution  of  2.'),  50,  25  and  an  overall  fuel  fuel-air  ratio  of  0.050.  the  local  fuel-air 
ratios,  top  to  bottom,  become  0.032.  0.113  and  0.032.  The  central  element  is  thus  operating  in  a 
very  rich  mode.  This  wide  distribution  requires  that  analysis  of  the  overall  combustion  efficiency 
must  be  performed  as  a multiple  streamtube  analysis. 

Utilizing  this  approach,  the  400‘’F  inlet  temperature  data  from  Run  9.01  were  analyzed. 
Figure  57  presents  the  results  of  this  study.  The  predicted  results  are  in  good  agreement  with  the 
test  data  with  regard  to  the  peak  efficiency  fuel-air  ratio  and  the  trend  of  efficiency  versus  f/a 
The  predicted  values  are  lower  than  the  test  data,  however.  The  major  problem  in  making  these 
predictions  b the  low  value  of  the  flame  speeds  at  the  fuel-air  ratio  extremes  generated  by  the 
fuel-air  distribution  problem.  As  a result,  the  results  are  quite  sensitive  to  the  value  of  laminar 
flame  speed. 

! 

This  value  is  currently  read  from  data  similar  to  Figure  58.  For  a typical  data  point  of  Run 
9.01,  the  fuel-air  ratios  are  calculated  to  be: 


Zone 

f/a 

Equivalence  Ratio 

Overall 

0.060 

0.882 

Top  F/H 

0.037 

0.544 

Center  F/H 

0.140 

2.059 

Lower  F/H 

0.039 

0.573 

From  Figure  58,  the  values  of  flame  speed  for  the  outer  two  zones  are  at  the  lower  end  of  the 
available  data,  but  the  central  zone  value  requires  an  educated  guess.  Since  this  one  zone 
comprises  22%  of  the  airflow  in  the  test  rig,  the  error  produced  by  this  one  zone  may  swing  the 
results  by  easily  10%  tie- 
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For  some  other  data  points  the  reverse  situation  is  true,  and  the  outer  two  zones  have  ver>’ 
lean  values,  and  the  local  equivalence  ratios  are  in  the  0..t  to  0.4  range.  Again,  the  flame  speed 
values  require  estimates.  For  example: 

Zone  Local  f/a  Elqui valence  Ratio 


Overall 

0.041 

0.602 

Top  F/H 

0.025 

0.367 

Center  F/H 

0.105 

1.544 

Lower  F/H 

0.021 

0.313 

The  results  of  Figure  57  were  generated  by  a three  streamtuhe  analysis  of  the  data, 
extending  the  flame  speed  cur\e  down  to  0.2  equivalence  ratio  and  up  to  2.0.  A symmetric 
behavior  of  the  flame  speed  versus  fuel  air  ratio  curve  about  a 1.1  equivalence  ratio  was  used. 

A similar  result  is  shown  in  Figure  59  for  the  data  from  Run  7.01  at  400'’F  inlet  and  120  to 
140  fps  approach  velocity.  The  comments  drawn  are  similar  to  those  from  Run  9.01. 

This  section  has  demonstrated  the  agreement  between  the  major  tenets  of  the  model  and  the 
experimental  results  of  th®  test  program.  Good  quantitative  agreement  was  obser\’ed  for 
predicted  and  observed  combustion  efficiency  data,  including  the  capability  to  predict  the  results 
of  fuel  distribution  changes.  This  latter  effect  is  observed  by  the  agreement  between  the  re.sults 
of  Run  7.01  and  the  predictions. 

b.  Specific  Results 

The  following  sections  will  address  the  specific  evaluations  of  geometric  and  aero- 
thermodynamic  variations  which  were  performed  during  the  test  program.  The  effects  of  fuel 
distribution  end  r verall  fuel-air  ratio  have  been  discussed.  As  in  the  previous  section,  this  report 
will  be  concerned  with  the  effects  on  overall  combustion  efficiency.  The  analysis  of  the  impact  on 
rumble  amplitudes  is  presented  in  the  report  on  the  Lo-Frequency  Augmentor  Instability  Study, 
Contract  F33615-76-C-2024. 

(1)  Heat  Addition  to  the  Wake 

The  wake  reaction  model  predicts  an  increase  in  wake  efficiency  and  a decrease  in  rumble 
gain  factor  if  heat  is  artificially  added  to  the  wake.  This  is  caused  by  an  increase  in  the  effective 
kinetic  reaction  rate.  On  the  test  rig  this  was  provided  by  an  oxyacetylene  torch  which  was 
mounted  to  feed  the  wake  of  the  central  flameholder. 

A problem  arose  during  this  test  in  that  a limitation  on  the  available  torch  flowrate 
prevented  achieving  the  levels  of  heat  addition  which  were  desired.  The  torch  heat  rate  of 
0.8  Btu/sec  to  the  wake  is  equivalent  to  a 15°F  rise  in  effective  inlet  temperature  or  a 3''r  increase 
in  reactor  efficiency  at  a 400°F  nominal  inlet  temperature.  At  this  level  of  heat  addition,  the  data 
scatter  normal  to  efficiency  measurements  exceeds  the  level  of  anticipated  improvements. 
However,  several  data  points  were  run  during  the  baseline  tests  with  the  52''r  blockage  and  the 
35*^^  blockage  flameholders  where  back-to-back  effects  were  noted  at  zero  and  full  torch  flowrate. 
These  results  are  shown  in  Table  4.  Figure  60  shows  the  results  over  a fuel-air  ratio  excursion  for 
35'^r  blockage  (Run  11.02).  In  general,  the  torch  resulted  in  a small  but  positive  gain  in  efficiency. 

This  result  is  significant  in  that,  although  small,  it  is  consistent  and  correct  in  direction  and 
magnitude  with  the  projected  effect.  Unfortunately,  the  direct  impact  on  rumble  amplitude  could 
not  be  determined  with  this  slight  heat  addition. 
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Average  Combustion  Efficiency 


TAitl.K  4 

KFFKrr  OF  WAKK  HKAT  ADDITION  ON  KFFICIKNCY 


Inlrt 

Temp  C^F) 
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Figure  60.  Effect  of  Heat  Addition  on  Overall  ('o<nbu‘::r  n. 

Efficiency  Run  11.02.  350  Bk>ckag‘'  Flamer.'  !der 


(2)  Effect  of  Available  Combustion  Length 

The  effect  of  the  available  axial  length  after  the  flameholders  was  examined  by  two  runs. 
Run  9.01  had  about  80  in.  available  for  post-flameholder  reaction  and  Run  10.01  reduced  this  to 
about  48  in.  These  teats  were  run  with  a 52''r  blockage  flameholder  and  a 40-in.  fuel  mixing  length 
at  200  and  400°F  inlet  temperatures.  During  this  test  .sequence,  the  fuel  was  observed  to  be 
trapped  by  the  volume  of  the  ports  to  the  observation  windows,  and  flameholding  occurred  at  that 
location  as  well  as  the  normal  flameholder. 

Additionally,  the  shortened  available  combustion  zone  ler<gth  pr  oduced  severe  sensitivity, 
and  rumble  or  unsteady  combustion  was  prevalent  over  lae  majority  of  the  test  sequence.  As  a 
result  of  these  problems,  no  useful  information  could  really  be  obtained  with  regard  to  efficiency. 

The  combustion  model  predicts  that  a shorter  length  results  in  a steeper  axial  temperature 
gradient  near  the  nozzle  inlet  and  lower  overall  efficiency.  The  steeper  gradient  should  have 
re.sulted  in  an  increased  sensitivity  of  the  average  inlet  temperature  to  pre.s.sure  and  velocity 
oscillations.  The  gross  system  behavior  certainly  agreed  with  this  qualitative  prediction. 

(3)  Effect  of  Inlet  Conditions 

The  effects  of  inlet  static  pressure,  temperature  and  approach  velocity  are  shown  in  Figures 
61  through  63.  The  effect  of  inlet  temperature.  Figure  61,  is  essentially  in  agreement  with 
predictions.  The  degree  of  increase  in  efficiency  is  lower  than  was  expected,  partially  due  to  the 
manner  of  the  tests  where  a fixed  level  of  inlet  flow  Mach  number  was  set.  The  fixed  Mach 
number  resulted  in  a I4^r  increase  in  absolute  velocity  at  a fixed  Mach  number  between  200  and 
400'’F  which  reduces  the  gain  due  to  the  increased  reaction  rates  at  the  elevated  temperature.  The 
remainder  is  unexplained  at  present. 

Although  the  data  is  sparse  for  a fuel-air  ratio  effect,  there  did  not  appear  to  be  any 
significant  effects  due  to  the  400  versus  200'’F  inlet  temperature.  Analysis  of  the  variations  in  the 
droplet  and  collected  film  vaporization  rates  predicted  only  about  a 10%  shift  in  vaporization 
rate. 


The  inlet  velocity  effect,  Figure  62,  is  also  in  the  anticipated  direction.  The  increase  in 
velocity  reduces  the  effective  flameholder  wake  reaction  efficiency  and  the  effective  flame 
spreading  angle.  These  effects  are  modeled  in  this  manner  in  the  analysis,  and  the  results  appear 
to  be  of  the  correct  magnitude. 

The  pressure  effect.  Figures  63  and  64,  is  interesting.  Although  the  data  are  quite  scattered, 
the  general  trend  is  a reduction  of  several  percent  in  overall  efficiency  when  the  pressure  is  raised 
from  15  to  21  psia.  This  trend  is  reversed  over  several  data  points  of  Run  8.01  in  Figure  64. 

The  predicted  effect  of  pressure  variation  on  the  fan  duct  combustion  process  has  offsetting 
parts.  The  vaporization  potential  at  fixed  liquid  temperature  is  reduced  as  the  static  pressure  is 
increased,  but  the  convective  heat  flux  increases  which  raises  the  liquid  temperature.  This 
increase  in  liquid  temperature  partially  or  completely  offsets  the  increased  pressure  effect  of 
reduced  vaporization  potential  for  these  test  data.  The  predicted  effect  is  a slight  increase  in  the 
overall  efficiency.  The  data  show  a slight  decrease,  which  may  be  scatter,  but  is  not  understood. 
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Figure  63.  Effect  of  Inlet  Static  Pressure  on  Efficiency 
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Figure  64.  Effect  of  Static  Pressure  on  Efficiency 
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(4)  Effect  of  Flameholder  Width  and  Blockage 

The  structure  of  the  test  apparatus  did  not  allow  an  independent  variation  in  flameholder 
width  and  blockage,  i.e.,  there  is  a fixed  total  flow  area  and  number  of  flameholders.  This 
coupling  of  width  and  blockage  produces  offsetting  effects.  The  increased  width  produces  an 
increased  wake  volume  but  also  increases  the  recirculation  rate.  The  increased  blockage  ratio 
generates  a higher  turbulence  level  at  the  flameholder  lip  and,  hence,  increased  flame  speed,  but 
also  increases  the  flameholder  lip  velocity  which  reduces  the  transverse  flame  spreading  rate. 

The  predicted  overall  efficiency  as  a function  of  width  and  blockage  was  shown  in  Figure  34. 
This  showed  a decrease  in  efficiency  of  about  lO^r  for  an  increase  from  ST.Sff  blockage  to  SO'V . 
The  test  data.  Figure  65,  shows  about  lOCf  reduction  in  i;c  for  3.5Cr  blockage  to  52fr  blockage  ratio. 
Although  the  amount  measured  is  slightly  lower  than  predicted,  the  agreement  is  felt  to  ue 
excellent.  This  section  of  the  model  is  fel  '.  to  be  in  an  accurate  form. 

(5)  Effect  of  Flameholder  Apex  Angle 

The  geometric  variation  of  flameholder  apex  angle  was  evaluated  at  the  3.5rr  blockage  level. 
The  standard  apex  angle  design  is  1.20-in.  wide  and  1.20-in.  for  a 53  deg  apex  angle.  The 
increased  angle  design  was  1.20-in.  wide  and  0.60-in.  deep  for  a 90  deg  included  angle.  This 
manner  of  altering  the  apex  angle  holds  width  and  blockage  con.stant  but  reduces  the  flameholder 
surface  area  by  This  change  will  thus  alter  the  rate  of  .surface  vaporization.  The  increased 
apex  angle  also  slightly  increases  the  rate  of  fuel  collection  onto  the  surface  and  the  rate  of 
gaseous  recirculation  into  the  wake  and  the  wake  volume. 

The  test  results  are  shown  in  Figure  66  for  400'’F  inlet  temperature  and  0.07  approach  flow 
Mach  number.  The  results  are  very  interesting  in  that  the  90  angle  flameholder  efficiencies  do  not 
exhibit  the  peak  welt  below  stoichiometric  that  was  ob.served  with  the  53  angle  design.  AKso,  the 
overall  levels  of  efficiency  appear  lower  than  with  the  smaller  apex  angle.  Comparison  of  Run 
11.02  and  12.01  at  22  psia  static  pressure  shows  that  the  peak  efficiencies  are  very  chwe  although 
Run  11  M nme  cl-city  oliould  have  been  tower. 

Analysis  of  the  effect  of  the  reduced  available  surface  area  produces  about  40%  surface 
vaporization  for  the  53  deg  angle  and  25%  for  the  90  deg  angle.  The  air  recirculation  constant,  K„ 
increases  from  23  to  25%,  and  the  collection  constant  increases  from  85  to  87%.  The  resultant 
effect  on  the  wake  fuel-air  ratio,  for  32%  preflameholder  vaporization  becomes: 

-|-= -Kl-ft)-^  (83) 

53  deg  angle:  (0.32)  + (0.68) = 1.325  (84) 

90  deg  angle:  (0.32)  + (0-68) = 0-912  (85) 

The  predicted  approach  fuel-air  ratios  for  peak  wake  reaction  are  thus: 

53  deg  angle:  f/a  = 0.051 

90  deg  angle:  f/a  = 0.076 

In  Figure  66  this  shift  was  definitely  observed. 
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igute  65.  Effect  of  Flameholder  Width  and  Blockage  Ratio  on  Combustion  Efficiency 


Overall  Fuel-Air  Ratio 


Figure  66.  Effect  of  Apex  Angle  on  Efficiency 

The  reduction  in  peak  efficiency  level  ia  not  in  agreement  with  the  predictions,  which 
expected  a small  gain.  This  error  ia  moat  likely  due  to  the  effect  on  the  larger  amount  of 
nonvaporized  fuel  which  leaves  the  flaineholder.  i 

The  static  pressure  influence  may  be  observed  agaii:  as  in  earlier  comments.: 

(6)  Effect  of  Free-Stream  Turbulence  Level 

The  impact  of  the  turbulence  level  on  the  efficiency  level  was  evaluated  by  the  u.s»>  of 
turbulence  generating  screens  placed  upstream  of  the  fuel  injector  spraybars.  This  screen 
generated  turbulence  of  a small  scale  suitable  for  interaction  with  the  recirculation  zone  shear 
layers  and  the  downstream  flame  propagation  ra'e. 

The  increased  turbulence  level  produces  two  basic  results  within  the  framework  of  the 
current  model.  The  increased  turbulence  level  produces  a higher  effective  recirculation  zone  shear 
layer  mass  exchange  rate.  This  increases  the  effective  unit  aerodynamic  loading  and  should 
reduce  the  wake  reaction  efficiency.  The  turbulence  also  increases  the  effective  turbulent  flame 
propagation  rate  downstream  of  the  flameholder. 

The  test  results  are  shown  in  Figure  67  for  400'’F  inlet  temperature  and  20-22  psia  static 
pressure.  The  increased  turbulence  level  produced  an  apparent  increase  of  7 to  10%  in  overall 
efficiency.  Apparently,  the  effect  of  increased  flame  spreading  rate  is  dominant  over  the  wake 
efficiency  reduction  in  controlling  overall  efficiency. 

Earlier  analysis  of  the  impact  'f  free-stream  turbulence  was  presented  in  the  model 
development  section.  The  trends  are  in  general  agreement  with  the  test  results.  For  example,  in 
Figure  35,  at.  increase  of  20%  in  turbulence  i.e.,  4 to  4.8%,  resulted  in  an  increase  in  efficiency 
from  70  to  76%.  The  test  data  at  0.055  overall  fuel-air  resulted  in  an  efficiency  increase  from  70 
to  78%.  These  cases  were  selected  to  yield  the  same  base  efficiency,  and  only  the  differences  are 
meaningful. 


Figure  67.  Effect  of  Turbulence  Level  on  Efficiency 

The  analysis  assumed  that  no  decrease  in  wake  efficiency  occurred.  In  reality  the  wake 
decreases  due  to  the  higher  mass  recirculation  rate  wh'ch  should  result  in  a lower  efficiency 
increase.  However,  the  screens  also  produced  a turbulence  increase  greater  than  the  20‘'f-  relative 
increase  used  in  this  analysis.  Thus,  the  match  between  data  and  predictions  is  quite  good. 

(7)  Effect  of  Inclined  (Drafted)  F ameholders 

The  effect  of  inclined  flam -holders  was  evaluated  on  Run  13.01  with  12-in.  wide  gutters 
which  were  angled  aft  at  a 22 ' i Jeg  angle  relative  to  the  spraybars.  Due  to  the  construction  of  the 
test  apparatus,  the  three  fiameholders  were  not  parallel  but  rather  staggered,  since  tbe 
attachment  port  for  the  central  flameholder  was  on  one  side  of  the  ducc  while  the  upper  and  lower 
flameholder  ports  were  on  the  opposite  side.  The  stagger  is  shown  in  Figure  68. 

The  flameholder  array  had  a projected  blockage  of  36%  and  1.2-in.  gutter  width.  The  swept 
ends  were  displaced  5 in.  downstream  resulting  in  13-in.  maximum  fuel  preparation  distance. 

Several  theories  have  been  advanced  for  why  drafted  gutters  should  or  should  not  improve 
performance.  Two  major  theories  are: 

1.  Sweeping  the  tips  of  the  flameholder  ters  downstream,  without  moving 
the  fuel  injector  spraybars,  increased  he  physical  separation  distance 
between  injection  and  the  flameholder.  1 his  will  result  in  a higher  level  of 
fuel  vaporization  and,  thus,  less  dependence  on  ♦he  surface  vaporization. 

Based  on  the  analysis  of  the  controlling  processes,  reduced  surtace  vapor- 
ization should  reduce  the  wake  fuel-air  ratio  and  extend  the  rich  limit.  For 
the  change  in  distance  as  tested  here,  the  difference  is  small. 

2.  The  drafting  causes  reduced  aerodynamic  blockage  which  generates  reduced 
turbulence  and  reduced  wake  aerodynamic  loading.  The  previous  test  results 
for  reduced  geometric  blockage  showed  this  to  result  in  increased  combustion 
efficiency. 
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Figure  68.  Inclined  Flcmeholders  vs  Normal  Location 

The  lest  results,  Figure  69,  show  a gain  of  5 to  in  overall  combustion  eTiciency.  however, 
the  comparison  data  from  Run  11.02,  also  35%  blockage,  are  at  a slightly  higher  approach  velocity 
which  should  reduce  its  efficiency  slightly.  From  these  data  the  conclusion  is  that  drafting,  by 
itself,  has  little  effect  on  efficiency. 

A three  point  test  was  run  at  T = 400®F  and  15  psia  to  evaluate  the  efficiency  fall-off  with 
approach  velocity.  The  test  results  are  shown  in  Figure  70.  The  severe  fall  in  efficiency  with 
velocity  is  as  the  analysis  predicts.  The  equivalent  data  from  the  undrafted  gutters  are  shown  for 
comparison.  No  discernible  effect  is  noted  which  suggests  that  the  wake  efficiency  versus  flame 
speed  trade  is  equal.  Unfortunately,  a velocity  blowout  limit  was  not  run  to  isolate  the  influences. 
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c.  Csncluslons 


The  agreement  with  the  experimental  results  which  the  two-phase  fuel  combustion  mode! 
has  generated  is  ver>’  encouraging.  The  major  conclusions  from  the  experimental  program  are  as 
follows; 

• The  major  thesis  of  the  two-phase  fuel  model  has  been  verified. 

• External  piloting  is  required  to  maintain  combustion  at  flow  Mach  numbers 
beyond  0.12  at  200  to  400“F  with  liquid  JP  fuels.  The  turbofan  engine  does 
this  by  transferring  core  engine  hot  gases  to  the  wake  of  the  fan  stream 
gutters. 

• Fuel  vaporization  processes  have  a major  impact  on  the  combustion  process. 

• Rumble  is  not  uniquely  tied  to  the  flameholder  near-wake  but  is  influenced 
by  it  insofar  as  the  near-wake  influences  the  flame  spreading  rate. 

• Rumble  tendencies  follow  the  classic  lines  of  combustion  instability  behavior. 

That  is,  any  change  which  reduces  the  velocity  or  pressure  sensitivity  of  the 
combustion  process  or  increases  the  system  damping  factor  will  reduce  the 
instability. 

• The  largest  rumble  producing  influence  which  is  not  predicted  by  the 
combustion  portion  of  the  model  is  the  .separation  distance  effect  of 
flameholder  to  fuel  injector.  Ihis  influence  is  a result  of  the  dvTiam.lcs  of 
liquid  fuel  oscillations  versus  air  oscillations  and  the  time  delay  between 
injection  and  stabilization. 

• The  multistreamtube  approach  to  modeling  fuel  distribution  effects  produces 
reasonable  results.  The  analysis  of  the  test  rig  ba.sed  on  three  streamtubes 
agreed  fairly  well. 

d.  Experimental  Mechanisms 

The  specific  conclusions  regarding  rumble  drivers  are: 

e Rumble  is  an  oscillation  in  the  overall  heat  relea.se  rate  of  the  augmentor 
produced  by  variations  in  the  flame  propagation  rate. 

e The  downstream  flame  propagation  rate  is  dependent  on  the  flameholder 
wake  reaction  efficiency.  Changes  which  alt»»r  this  efficiency  produce 
alterations  in  overall  heat  relea.se  rate,  and  sensitivity  of  the  proce.ss  to 
velocity  or  pressure  oscillations. 

• A fuel-air  ratio  distribution  other  than  uniform  is  destabilizing. 

• Any  change  which  reduces  the  flame  spreading  rate  is  generally  destabilizing: 

Increased  velocity 

Reduced  temperature 

Reduced  wake  size  (residence  time) 

Rich  wake  fuel-air  ratio. 
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• The  following  postulated  mechanisms  were  rejected: 

— Fuel  injector  flowrate  oscillations 
— Wake  vortex  shedding 
— Fuel-air  ratio  oscillations. 

4.  PHASE  II  — MODEL  EVALUATION 

a.  Task  I — Comparison  of  Model 

The  efforts  under  this  task  were  essentially  grouped  into  two  categories: 

• Computerization  of  Model 

• Comparison  with  other  Models  and  Data. 

The  two  tasks  were  performed  simultaneously  as  the  computer  analysis  was  refined  and 
predictions  compared  with  available  data. 

(1)  Computerization  of  the  Model 

The  analyses  performed  during  Phase  I were  assembled  into  a cohesive  computer  program 
which  is  capable  of  analyzing  either  the  core  or  fan  duct  stream  of  a conventional  mixed  tlow 
turbofan  augmentor.  The  program  requires  as  input,  a description  of  the  geometry  of  the 
augmentor  stream  and  a description  of  the  physical  operating  conditions  to  be  analyzed. 

The  operational  logic  is  such  that  each  case  run  represents  a single  geometry  and  opera'ing 
point.  The  program  will  generate  a curve  of  the  axial  average  combustion  efficiency  and  final 
augmentor  section  performance. 

The  program  operational  logic  is  shown  in  Figure  71.  and  the  subroutines  are  defined  in 
Table  5. 

In  addition  to  the  augmei.tor  duct  axial  efficiency  curve,  the  program  evaluates  the  wake 
conditions  for  fan  stream  flameholders.  Based  on  the  technical  development  previously 
discussed,  the  duct  flameholder  wake  composition  is  defined  as: 


S,  + (1-d.)  • 


Ml 

K, 


(80) 


The  computer  code  evaluates  these  coefficients  as  well  as  the  degree  of  wake  reaction 
efficiency,  i;,,  at  the  wake  fuel-oil  ratio,  0,  and  the  wake  aerodynamic  loading.  A/v„,,b 


The  solution  procedure  for  the  droplet  vaporization  and  collection  onto  the  flameholder 
proceeds  directly.  The  recirculation  rate  solution  is  also  accomplished  easily.  A probler.n  arose 
with  the  solution  for  the  surface  vaporization  rate,  and  the  wake  reaction  efficiency.  ?(..  The 
problem  was  du’  to  the  interdependent  natur-'  of  and  i;„.  The  vaporization  depends  strongly 
on  the  wake  temperature  which  is  a function  of  the  wake  fuel-air  ratio  and  aerodynamic  loading. 
Since  the  wake  fuel-air  ratio  is  a function  of  0,.  a direct  explicit  .solution  is  not  possible. 


Attempts  to  write  an  explicit  algorithm  for  the  solution  did  not  produce  the  desired  rcsu'*s. 
This  failure  was  primarily  due  to  the  nature  of  the  wake  kinetics  equations  which  are  also  implicit 
functions  of  the  aerodynamic  loading  in  terms  of  wake  efficiency  and  fuel-air  ratio,  e.g. 


fen  (0,  i/w) 
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Col/ect 


Figure  71.  Basic  Computer  Program  Logic  Diagram 
96 


fD  I40(rs 


TABLE  5 

PROGRAM  SUBROUTINE  DESCRIPTION 


Subroutine 

Major  Operations 

Inject 

Seta  up  droplet  size  groups 

Injection  vaporization  level 

Accel 

Droplet  vaporization  and  acceleration  from 
sprayring  to  flamebolder 

Collect 

Droplet  collection  onto  flamebolder  surface 

Recirc 

Wake  recirculation  rate  and  aerodynamic  load- 
ing 

Beta  3 

Liquid  vaporization  rate  from  surface  film 

Wake 

Wake  kinetics  solution 

Flame 

Two-phase  fuel  flame  .speed 

Flame  2 

i 

Gaseous  flame  speed 

The  resultant  equation  set  which  must  be  solved  is  thus: 

• i 

« = + 9 (1-/9.) 

M. 

Ki 

0t  = fen  (T„  di.  9) 

1 

i 

Tw  = T,  + ATid„ 

; 

AT,*.,,  = fen  (T„  0) 

(88) 

(89) 

(90) 

(91) 


= fen  (T«,  0,  A/V„p*)  (92) 

Where  the  i;,  function  has  the  form  shown  earlier  in  equation  87. 

The  known  quantities  are  9,  j3,,  K,,  Ta,  AA^„p*,  and  the  ATw,.,  function.  The 

solution  should  yield  T,,  <p,  rjy,.  The  required  flexibility  of  the  analysis  should  also  be  able  to 
recognize  the  possibility  that  a solution  does  not  exi.st.  For  example,  over  the  range  of  wake  fuel- 
air  ratios,  a maximum  value  of  kinetic  capability  exists.  The  analysis  must  recognize  an  input  set 
which  produces  an  aerodynamic  loading  in  excess  of  this  value  as  an  augmentor  blowout. 


Another  failure  might  occur  when  the  compositional  solution  (Equation  88)  does  not  yield 
a w.i!:e  fuel-air  ratio  which  is  possible  from  the  kinetics  solution  (Equation  87).  Typically,  this 
occurs  as  a rich  blowout.  Both  of  these  failures  are  a direct  con.sequence  of  the  input  values.  The 
program  is  written  to  analyze  a situation  which  exceeds  the  physical  limits  for  that  particular 
geometry  and  flow  field. 
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(2)  Computer  Program  Description 


The  computer  code  is  currently  set  to  analyze  a conventional  turbofan  augmentor  with  vee- 
gutter  flamcholders.  The  generalized  geometry  which  is  analyzed  is  shown  in  Figure  73.  The 
liiiuid  fuel  is  injected  through  a spraybar  located  a specified  distance  upstream  of  the 
flameholder.  The  combustion  process  is  stabilized  by  the  bluff  body  recirculation  zone,  and  a 
turbulent  flame  sheet  propagates  into  the  approach  fuel-air  mixture.  The  geometry  is  two- 
dimensional  and  would  repre.sent  one  .streamtube  of  a multiple  stream  augmentor  system. 

The  various  analyses  which  were  developed  previously  and  referenced  in  the  following 
paragraphs  are  written  as  subroutines  in  the  computer  program.  This  results  in  a modular 
program  with  a supervi.sory  MAIN  executing  subroutines  as  required.  This  approach  allows  easy 
modification  of  the  various  analyses  without  di.sturbing  the  overall  program  operational  logic. 

The  program  first  reads  the  input  in  NAMELIST  format.  From  the  input  the  two- 
dimensional  model  format  is  set  up.  Total  air  and  fuel  flowrates  are  calculated  and  the  core 
stream  or  duct  stream  option  is  exercised.  The  following  description  details  the  duct  stream 
analysis  procedure. 

From  the  fuel  conditions  in  the  spraybar  and  the  flow  field  conditions  in  the  duct,  the  degree 
of  flash  vaporization  which  occurs  during  injection  is  calculated.  This  percentage  of  fuel  is 
allocated  to  the  initial  vapor  phase. 

The  fuel  which  remains  liquid  is  assigned  to  five  equal  ma.ss  flowrate  groups.  These  groups 
are  each  assigned  a mean  droplet  diameter.  The  droplet  diameters  generated  by  the  program 
repre.sent  the  spray  formation  characteristics  of  a variable  area  pintle  spraybar.  The  controlling 
size  parameter  is  the  fuel  injection  pressure  drop.  These  data  are  empirical  from  Pratt  & Whitney 
Aircraft  data.  These  five  groups  thus  repre.sent  the  flowrate  versus  size  distribution  unique  to  this 
spraybar  type.  If  a different  type  of  injector  is  to  be  analyzed,  the  droplet  sizing  subroutine  must 
be  rewritten. 

The  program  then  performs  a finite  difference  solution  to  the  droplet  acceleration  and 
vaporization  equations  by  selecting  a small  time  step  and  solving  for  the  deltas  of  liquid  velocity, 
temperature  and  a delta  vaporized  from  the  drag  and  vaporization/heating  equations.  The  axial 
travel  of  the  liquid  droplet  during  this  time  increment  is  calculated  from  the  initial  and  final 
liquid  velocities.  The  analysis  continues  this  small  time  step  solution  until  the  axial  distance 
value  equals  the  spraybar  to  flameholder  separation  distance.  The  acceleration/vaporization 
mfxlel  assumes  that  the  isolated  droplet  vaporization  rate  exceeds  the  rate  in  spray  clouds  by  a 
factor  of  two.  The  calculation  sums  the  amount  of  fuel  which  vaporizes  from  all  the  droplet  size 
groups  and  adds  this  to  the  amount  of  fuel  vaporized  during  injection.  The  sum  of  these 
represents  the  total  term  in  the  wake  compositional  equation. 

The  analysis  next  calculates  the  percentage  of  each  of  these  size  groups  which  is  collected 
by  impingemc-iit  onto  the  flameholder  surface.  The  collection  rate  is  calculated  for  each  size  group 
individually  utilizing  the  droplet  diameter  which  exists  for  each  group  after  the  vaporization 
calculations.  The  total  collected  liquid  flowrate  is  summed  from  the  collection  efficiency  for  each 
group  and  the  liquid  flowrate  which  exists  in  each  group  after  the  vaporization  calculations.  The 
collection  rate  calculation  for  each  droplet  size  group  is  evaluated  from  droplet  trajectory 
an.dy.ses  which  were  performed  and  correlated  against  the  geometric  variables  of  the  flameholder 
system.  The  details  of  this  were  presented  in  the  earlier  Phase  I results. 
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The  wake  recirculation  rate  is  next  evaluated,  again  following  the  earlier  reported  analyses. 
The  influence  curves  were  reduced  to  a series  of  equations  or  curve  reading  routines  within  the 
RECIRC  subroutine.  This  subroutine  evaluates  the  aerodynamic  loading  of  the  idealised 
recirculation  zone  based  on  the  empirical  data  earlier  pre.sented.  The  recirculation  coefficient,  K,, 
is  calculated,  and  the  loading  is  stored  for  tran.sfer  to  the  wake  kinetics  subroutine. 

At  this  point  the  solution  proceeds  along  the  parallel  paths  of  exercising  the  flamcholdcr 
surface  vaporization  subroutine,  BETA  3,  and  the  wake  kinetics  subroutine. 

For  the  solution  of  BETA  3,  values  of  wake  temperature  are  assumed  for  every  200°F 
increment  from  1000  to  5000“F.  For  each  of  these  values,  the  heat  flux  from  the  recirculation  zone 
through  the  flameholder  into  the  liquid  film  is  calculated.  The  surface  vaporization  is  calculated 
utilizing  a 20-step  finite  difference  solution  to  the  forced  convection  vaporization  problem.  The 
solution  technique  is  the  same  as  presented  earlier.  A 20-step  solution  was  found  nece.ssary  for 
convergence  of  numerical  accuracy  and  to  assure  a smooth  evaluation  as  the  vaporization  rate 
passes  through  to  lOO^r.  A check  is  made  in  the  calculation  for  the  ratio  of  heat  flux  to  the  liquid 
latent  heat  at  collection  conditions.  Whenever  the  ratio  exceeds  unity,  the  vaporization  is  set  at 
lOOfi. 

For  each  value  of  wake  temperature,  the  calculated  vaporization  rate  is  used  in  conjunction 
with  the  previous  compositional  coefficients  to  define  a wake  vapor  phase  fuel-air  ratio.  This 
temperature  versus  fuel-air  ratio  array  is  stored  for  later  use. 

The  program  transfers  next  to  the  WAKE  subroutine  for  the  kinetics  solution  utilizing  the 
recirculation  rate  from  RECIRC  and  the  Longwell  reactor  model  pre.sented  earlier.  The 
subroutine  is  exercised  over  a range  of  fuel-air  ratios.  The  fuel-air  ratios  are  started  at  0.020  and 
increased  by  0.0045.  The  solution  is  stored  in  an  array  of  efficiency  and  wake  temperature  versus 
wake  fuel-air  ratio  between  the  lean  and  rich  limits. 

These  limits  are  evaluated  at  the  level  of  aerodynamic  loading  for  the  case  being  analyzed. 
The  lean  limit  is  the  first  fuel-air  ratio  where  a kinetic  solution  is  found.  The  rich  limit  is  found 
by  a fine  grid  search  backwards  from  the  first  fuel-air  ratio  which  fails  to  produce  a solution  on 
the  rich  side.  The  increased  accuracy  of  definition  of  the  rich  limit  was  found  necessary  since  the 
rich  blowout  is  the  significant  failure  mode  of  duct  flameholders. 

Once  this  array  is  generated,  the  program  searches  for  array  intersection  between  the 
BETA  3 and  WAKE  results.  If  none  is  found,  appropriate  failure  messages  are  printed  and 
execution  stops.  If  intersection  is  found,  the  convergent  re.sults  of  the  wake  compositional  solution 
are  printed. 

The  next  subroutine,  FLAME,  performs  the  finite  difference  solution  to  the  two- 
dimensional  flame  propagation  problem  presented  earlier.  The  flame  speed  base  value  is  altered 
by  the  level  of  v/ake  reaction  efficiency  and  by  the  percentage  of  liquid  fuel  which  is  sloughed  from 
the  flameholder  trailing  edge.  This  influence  accounts  for  two  real  effects  in  the  augmentor 
transferred  to  the  pseudo-two-dimensional  analysis.  The  wake  efficiency  is  assumed  uniform  and 
continuous  as  is  the  flame  sheet.  In  reality,  as  the  efficiency  decreases,  local  regions  are  generated 
where  flame  initiation  in  the  shear  layer  fails.  Decreasing  the  transverse  speed  of  the  idealized 
flame  sheet  is  a method  of  describing  this  effect  to  account  for  the  requirement  of  flame 
propagation  normal  to  the  duct  into  these  unignited  regions. 

Similarly,  the  sloughed  fuel  serves  to  quench  local  regions  of  ignition  in  the  shear  layers 
which  also  require  more  normal  flame  penetration.  The  inclusion  of  a vaporizntion/sloughing 
term  accounts  for  this  in  a two-dimensional  analysis. 
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Once  the  flame  is  initiated,  the  subroutine  performs  a finite  difference  solution  to  the 
transverse  propagation  into  the  free-stream  allowing  for  the  axial  variations  in  velocity, 
turbulence,  etc.  This  continues  until  the  exhaust  nozzle  is  reached. 

For  the  analysis  of  a core  stream  segment  the  procedure  is  simplified  greatly  due  to  the 
complete  fuel  vaporization.  The  program  sets  the  wake  fuel-air  ratio  equal  to  the  total  fuel-air 
ratio  and  performs  the  wake  kinetics  solution  at  this  value  and  the  calculated  value  of 
recirculation.  A gaseous  phase  turbulent  flame  is  initiated,  corrected  by  wake  efficiency,  and 
evaluated  downstream  as  before. 

(3)  Model  Comparison  With  Other  Models  and  Data 

The  model  for  the  stabilization  of  two-phase  fuel  at  low  inlet  temperatures  did  not  employ 
analog  data.  This  model  was  evaluated  against  the  available  literature  references.  The  analysis 
of  the  high  inlet  temperature  core  stream  stabilization  and  propagation  model  was  performed  by 
comparison  with  the  currently  used  turbojet  augmentor  efficiency  correlations  as  well  as  tbe 
available  body  of  gaseous  fuel  stabilization  literature. 

The  two-phase  stabilization  model  is  based  on  the  pioneering  research  efforts  of  Marion  in 
1952  (Reference  31).  In  his  studies,  Marion  generated  the  basic  concept  of  flame  stabilization 
under  conditions  of  flameholder  surface  vaporization.  The  testing  and  analysis  were  performed 
for  low  volatility  fuels  (primarily  diesel  fuel)  but  the  qualitative  behavior  of  the  situation  is 
comparable. 

The  major  consequence  of  reduced  vapor  pressure  fuels  on  the  process  of  surface 
vaporization  is  in  the  location  of  the  equilibrium  solution  to  the.  wake  vapor-phase  fuel-air  ratio. 
For  the  reduced  volatility  diesel  fuel  employed  in  Marion’s  studies  two  compositional 
intersections  were  calculated.  One  on  the  lean  side  was  rejected  as  a solution  due  to  its  lack  of 
stability.  The  solution  on  the  rich  side  was  stable,  however.  The  solutions  for  the  more  volatile  JP 
fuels  produce  a single  intersection  of  wake  formative  and  reaction.  This  solution  lies  on 
vaporization  coefficient  less  than  1.0.  The  value  may  be  on  the  lean  side,  if  the  surface 
vaporization  is  100%,  and  still  be  stable. 

Figure  72  is  a typical  compositional  result  generated  by  Marion.  Figure  74  is  a similar  result 
from  the  model  generated  here.  The  variation  in  the  intersection  behavior  is  the  result  of  the 
variations  in  surface  vaporization  rates  for  given  wake  heat  transfer  caused  by  fuel  volatility 
differences.  ^ 

Another  study  on  wake  flame  stabilization  with  liquid  fuels  was  performed  by  R.  M.  Gross 
(References  32  and  33).  In  this  study  the  liquid  fuel  was  injected  on  thejupstream  side  of  a semi- 
porous  cylindrical  flameholder.  Two  interesting  results  from  this  study  confirm  critical  sections 
of  the  analysis. 

First,  Gross’  experimental  program  confirmed  the  behavior  of  t’ne  surface  vaporization 
phenomenon  as  predicted  by  this  study.  The  experimental  efforts  employed  wake  compositional 
sampling  at  various  flowrates  of  the  surface  liquid  film.  The  observed  wake  behavior.  Figures  75 
and  76,  show  a plateau  in  composition.  This  was  observed  to  coincide  with  the  onset  of  liquid 
sloughing  from  the  surface  as  the  fuel  flowrate  was  increased.  Essentially  this  implies  that  the 
vaporization  rate  maximum  value  is  controlled  by  the  available  heat  flux  from  the  wake.  Any 
excess  fuel  is  merely  lost.  This  effect  is  not  linear  due  to  the  heat  transferred  to  the  sensible  heat 
of  the  lost  liquid.  In  fact,  t!;  3 data  may  be  read  as  exhibiting  a slight  recurvature. 
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Figure  74.  Flame  Stability  Criteria 
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The  effect,  as  predicted  by  the  rumble  combustion  model,  v.as  shown  in  Figure  19  as  the 
effect  of  fuel  collection  rate  on  the  surface  vaporization  rate.  The  decrease  in  percentage  as  the 
total  fuel  flowrate  is  increased  will  yield  the  same  behavior  as  seen  by  Gross.  The  quantitative 
shift  is  due  to  the  lower  air  velocities  used  by  Gross  wbicb  result  in  lower  mass  effect  rates  due 
to  the  reduced  mass  transfer  Nussett  number. 

The  second  major  point  of  the  model,  which  was  confirmed  by  Gross’  test  results,  is  the 
influence  of  the  approach  air  velocity  on  the  wake  composition  (Figure  77).  There  is  little  or  no 
influence  at  low  velocities  but  a very  rapid  increase  in  unreacted  composition  as  the  blowout  limit 
is  reached.  The  sharp  oxygen  concentration  increase  is  primarily  due  to  the  rapid  fall-off  in 
reaction  kinetic  efficiency.  The  comparison  results  (Figure  78)  show  a sharp  decrease  in  CO,  at 
the  same  time.  The  CO  concentration  shows  an  increase  at  first,  due  to  the  failure  of  the  relatively 
slow  CO  oxidation  reaction  and  a fall-off  near  blowout  due  to  the  failure  of  the  CO  formation 
kinetics.  This  data  tends  to  confirm  the  kinetic  model  of  the  wake  as  the  controlling  parameter 
as  did  Longwell’s  original  analysis  (Reference  13)  used  in  this  model. 

The  augmentor  model  was  exercised  against  development  experience  on  a mixed  flow 
turbofan  engine.  Predictions  wore  made  for  the  fan  duct  and  core  stream  efficiencies  as  a function 
of  altitude  at  a fixed  subsonic  flight  Mach  number.  The  influence  of  heat  addition  to  the  wakes 
of  the  fan  duct  flameholders  was  also  predicted,  since  this  was  studied  on  the  engine  and 
produced  significant  gains  in  the  rumble  onset  altitude. 


Figure  77.  The  Effect  of  the  Entrance  Air  Velocity  on  the  Eddy  Gas  Composition  for  the 
18  X 10~*  tbm/(in.)*  (min.)  JP-4  Wake  Flame:  H,0,  N,,  and  0,  Concentration  Profiles 
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Figure  78.  The  Effect  of  the  Entrance  Air  Velocity  on  the  Eddy  Gas  Composition  for  the 
18  X /0‘*  tbrn/dn.)  ( min.)JP’4  Wake  Flame:  COt,  CO,  and  H,  Concentration  Profiles 

The  results  for  the  core  stream  are  shown  in  Figure  79.  The  behavior  follows  the  classic  curve 
of  efficiency  versus  fiiel-air  ratio  for  gaseous  fuels.  The  peak  efficiency  fuel-air  ratio  does  not  occur 
at  stoichiometric  (0.068  for  JP-4)  but  rather  at  about  0.055.  This  shift  is  caused  by  the  vitiation 
of  the  turbine  exit  flow  due  to  main  burner  combustion.  For  the  cases  analyzed,  the  core  stream 
efficiency  is  quite  good  up  to  very  high  altitudes.  This  would  be  expected  at  this  level  of  approach 
temperature. 

The  analysis  for  the  fan  duct  stream  resulted  in  quite  different  behavior.  The  efficiency 
(Figure  80)  shows  a sharp  climb  from  the  lean  limit  and  a plateau  of  efficiency  for  some  range  of 
fuel-air  ratios  followed  by  a sharp  fall  towards  rich  blowout.  Also,  significant  decline  occurs  above 
60,000  feet  at  the  ana  zed  flight  Mach  number. 

The  peak  efficiency  fuel-air  ratio  is  approximately  0.035.  This  is  due  partly  to  the  relative 
wake  enrichment  from  the  two-phase  fuel  concept  and  partly  to  the  loss  of  some  combustion 
airflow  diverted  for  liner  cooling.  The  sharp  decline  past  the  0.050  fuel-air  ratio  is  due  to  the  wake 
nearing  its  rich  limit.  This  is  accompanied  by  a rapid  decline  in  reaction  capability  with 
increasing  fuel-air  ratio. 
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Figure  79.  FlOO  Core  Stream  Efficiency  Prediction 


Duct  Fuel-Air  Ratio 

n>  i«nM 

Figure  80.  FlOO  Fan  Duct  Stream  Efficiency 
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'I'wi)  factors  contribute  to  the  rapid  decline  in  efficiency  from  oD.diXl  to  nojMiO  h i t.  1 lu’ 
reduced  static  pressure  in  the  nu(,’mentor  increases  the  rate  ttf  fuel  vaporization  from  the  surhu  c 
of  the  flameholder  for  a Riven  level  of  vtake  heat  flux.  Also,  the  reduced  pressure  nsulis  m i 
decreti.se  in  the  vake  rich  limit  fuel-air  ratio,  'rhtis,  a richer  wake  with  ;i  le.ini'r  rich  limit  \ icl(l> 
earlier  blowcmts  at  higher  altitudes.  This  region  is  shown  in  more  dcttiil  in  Figure  81 . This  region 
is  the  same  altitude  region  where  the  developmetu  engine  experienced  rumhle  insttihiliiy  tit  fuel 
air  ratios  above  O.O.uO.  The  predictions  are  shown  below; 


Altitude 

Rumhle  Limit  f/n 

51,000 

o.orx) 

5;i,0()0 

0.058 

54,0(X) 

0.052 

5.0.UOO 

o.o:)8 

The  extreme  sensitivity  of  this  limit  underscores  a recurring  developmeni  engine  problem. 
During  evaluation  testing  for  rumble  limit  altitudes  or  fuel-air  ratios.  se\ere  data  scatter  is 
continuously  encountered.  The  fan  duct  behavior  is  such  that  small  shills  in  the  test  engine  trim 
or  fuel  schedule  within  normally  allowable  error  bands  will  generate  large  changes  in  the  limit 
points. 

The  predicted  behavior  of  the  engine  is  right  in  the  development  engine  diita  range  undi  r 
these  conditions  of  flight  Mach  number.  The  total  augmentor  efficiency  versus  fuel-air  ratio  at 
50,000  feet  altitude  is  shown  in  Figure  82.  These  predictions  agree  fairly  well  with  the  engine 
behavior  at  this  flight  point.  Again  the  severe  decrease  above  ahou*  0.050  fuel-air  ratio  is  evident . 
At  these  conditions,  the  general  experience  rumble  limit  was  in  the  0.052  range,  which  agues  well 
with  the  predictions. 

The  influence  of  wake  heat  addition  is  shown  in  Figure  8.'?.  A significant  increase  in  the  limi' 
fuel-air  •■atio  is  predicted  as  well  as  a slight  gain  in  overall  efficiency  at  lower  fuel-air  ratios.  'I'he 
dramatic  increase  in  limit  f/a  at  the  higher  altitudes  is  [primarily  due  to  the  increase  in  the 
allowable  rich  limit  wake  f/a.  The  heat  addition  greatly  increases  the  wake  kinetic  capability 
which  expands  the  rich  limit. 

Two  data  comparisons  arc  possible.  First,  the  data  from  the  experimental  program  run  for 
this  study  showed  an  increase  of  2 to  -Vi  in  augmentor  efficiency  with  direct  wake  heat  addition 
The  model  predicts  about  4 to  5*^’;  gain  for  the  engine  at  similar  conditions.  .Second,  the  engiiio 
development  experience  was  an  altitude  increase  of  .50(X)  to  10, WX)  feet,  at  0.055  f/a,  in  rumble 
limit  altitude.  The  model  predicts  appioximately  7000  feet.  This  heat  addition  on  th.e 
development  engine  resulted  in  problems  not  related  to  this  study  which  iirecluded  its  use. 

'I’he  basic  portions  of  the  model  were  previously  compared  to  the  test  data  from  the 
exirerimental  program  and  also  compared  to  a generalized  turbojet  engine  efficiency  corri'h'tion 
with  good  results.  'Fhese  comparisons  were  presented  in  the  earlier  section  of  this  re|)ort . Ty[)ical 
results  are  repeated  in  Figime  84. 
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Figure  81.  Fan  Duct  Efficiency  Near  the  Blowout  Limit 
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Figure  82.  FIOO  Augmentor  Predicted  Efficiency  vs  Fuel-Air  Ratio 


Figure  83.  Effect  of  Heat  Addition  on  Duct  Efficiency 
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b.  Task  II  — Recommend  Improved  Designs 

This  task  involved  two  related  elements  involving  extension  of  the  model’s  results  into  the 
area  of  augmentor  rumble  elimination.  These  elements  were: 

e Extend  the  model  to  other  concepts 

a Predict  the  stability  of  these  concepts. 

This  task  was  a relatively  minor  part  of  the  total  program  in  terms  of  expended  activity  but 
important  in  terms  of  identification  of  directions  for  improved  turbofan  augmentor  technology. 

The  model  formulated  was  derived  specifically  for  the  conventional  bluff  body  stabilized 
turbofan  augmentor  with  normal  turbulent  flame  propagation  controlling  the  ultimate  level  of 
augmentor  efficiency.  The  primary  cause  of  rumble  is  the  dynamic  response  of  these  processes  to 
fluctuations  in  either  inlet  conditions  or  fuel-air  ratio. 

This  similar  augmentor  response  logic  was  used  for  definirg  two  advanced  augmentation 
concepts.  The  analyses  for  these  concepts  are  less  rigorous  than  for  the  conventional  augmentor 
for  two  reasons; 

» Considerably  less  effort  was  expended 

• Less  is  known  about  the  controlling  processes  than  for  the  conventional  case. 

The  two  concepts  examined  here  are  the  swirling  flow  augmentor  and  the  VORPIX 
augmentor.  In  both  of  these  concepts  there  is  an  attempt  to  provide  a “hard”  pilot  rather  than  the 
“soft”  flameholder  wake  region  and  to  reduce  the  dependence  of  the  downstream  transverse  flame 
penetration  rate  on  the  fuel-air  ratio.  Both  of  these  changes  will  serve  to  reduce  the  augmentor 
efficiency  response  and  will  generate  increased  stability  limits. 

The  swirling  flow  augmentor  concept  generates  an  increa.sed  flame  speed  by  implementa- 
tion of  the  concept  of  bouyant  penetration  in  centrifugal  force  fields.  The  basic  concepts  are 
presented  in  Reference  34  and  subscale  rig  test  experience  in  Reference  35.  The  es.sential  elements 
of  an  augmentor  utilizing  swirling  flow  consist  of  inlet  flow  swirl  vanes,  an  outer  diamet  er  annular 
pilot  with  its  own  fuel  system  and  zoned  circumferential  fuel  sprayrings. 

The  outer  pilot  provides  a circumferentially  uniform  flame  initiation.  The  transverse  flame 
propagates  via  centrifugally  enhanced  turbulent  flame  speed  into  the  swirling  approach  flow.  The 
flame  speeds  achieved  by  this  mechanism  are  many  times  greater  than  from  normal  turbulent 
flame  penetrations,  fn  addition,  the  driving  potential  for  the  flamt  speed  is  the  bouyant  force 
caused  by  density  differences  between  the  burned  and  unburned  gas»s  rather  than  turbulent- 
kinetic  transport.  Thus,  the  variation  due  to  fuel-air  ratio  excursions  is  much  less. 

Figure  85  compares  the  predicted  variation  in  flame  speed  versus  fuel-air  ratio  for  a .swirling 
flow  and  a conventional  design.  The  curves  are  both  normalized  to  1.0  at  the  maximum  flame 
speed  fuel-air  ratio.  The  dramatic  difference  is  evident.  This  degree  of  change  will  reduce  the 
variational  influence  of  fuel-air  ratio  excursions  on  total  augmentor  efficiency. 

The  results  of  the  subscale  rig  testing  tended  to  verifv  the  improved  combustion  efficiency 
over  a wider  fuel-air  ratio  range.  It  also  underscored  the  very  strong  influence  of  the  fuel 
distribution  on  the  augmentor  performance.  Typical  results  are  shown  in  Figure  86.  This  test 
apparatus  utilized  four  zones  of  fuel  injection  with  zone  1 feeding  the  annular  pilot  and  zones  2, 

3 and  4 located  more  toward  the  centerline.  The  largest  efficiency  values  were  obtained  by 
maintaining  as  nearly  as  possible  locally  stoichiometric  conditions. 
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The  effect  of  inlet  velocity  during  testing  was  considerably  less  than  conventional 
flameholders  ns  shown  in  Figure  87.  The  variation  with  fuel-air  ratio  is  shown  in  Figures  88  and 
89  at  two  different  augnientor  lengths.  The  insensitive  nature  is  clearly  seen. 

The  anticipated  behavior  of  an  augmentor  design  utilizing  the  swirl  concept  would  be  an 
extremely  stable  configuration.  However,  there  are  several  (pialifying  factors;  the  major  ones  are 
delineated. 

• The  fuel  distribution  should  be  designed  to  produce  local  stoichiometry  as 
nearly  as  possible.  The  data  clearly  show  that  local  richness  will  produce 
efficiency  declines  severe  enougn  to  drive  rumble  instabilities 

® The  variational  effect  of  pressure  on  the  flame  speed  needs  experimental 
evaluation 

• The  pilot  must  be  sized  to  eliminate  any  significant  airflow  dynamic 
response.  If  the  pilot  responds,  the  effect  will  be  the  same  as  if  the 
flameholder  wake  is  forced  near  blowout 

• More  research  is  required  on  the  fundamental  processes  involved,  particular- 
ly the  effect  of  liquid  fuel  concentration  which  at  present  is  virtually 
unknown. 


Fif^urc  85.  Sirirlinf’  Flow  Flame  Speed  vs  Conventional  Flame  Speed 
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Open  Symbols  - L/D  = 1.414;  Closed  Symbols  - L/D  = 0.914 


O - Zone  2 Only 

□ - Zone  2 0 = 0.175;  Zone  3 0 Varied 

^ - Zone  2 0 = 0.339;  Zone  3 0 Varied 

^ - Zone  2 0 = 0.518;  Zone  30  = 0.134 

- Zone  2 0 = 0.328;  Zone  30  = 0.335;  Zone  4 0 Varied 

0 - Zone  2 0 = 0.397;  Zone  3 0 = 0.167;  Zone  40  = 0.339 

O - Zone  2 0 = 0.699;  Zone  4 0 Varied 

4 - Zone  2 0 = 0.463;  Zone  30  = 0.432 

A • Zone  2 0 = 0.486;  Zone  4 0 Varied 

^ - Zone  2 0 = 0.286;  Zone  40  = 0.168 


Figure  86.  Effect  of  Augmentor  LfD  on  Combustion  Efficiency 
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Combustion  Efficiency  Percent 


Augmentor  Equivalence  Ratio  x 10 

FD  14013f 

Figure  88.  Comparison  of  Gas  Sample  and  Choked  Nozzle  Calculated  Efficiencies 
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Figure  89.  Comparison  of  Gas  Sample  and  Choked  Nozzle  Calculated  Combustion  Efficiencies 

The  alternative  augmentor  design  to  which  the  rumble  model  concepts  were  applied  is  the 
P&WA  VORBIX  augmentor.  This  augmentor  concept  (Reference  36)  relies  on  a pilot  burner  fed 
by  high-temperature  engine  stream  airflow.  The  main  fuel  injection  is  into  the  pilot  discharge 
flow  where  vaporization  occurs.  This  fuel  rich  vapor  cloud  is  penetrated  later  by  vortices  of  the 
fan  duct  and  engine  airflows  which  are  driven  by  mechanical  swirling  devices.  This  produces  a 
central  spinning  core  of  cooler  air  surrounded  by  hot  fuel  rich  vapor. 

This  situation  produces  flow  instability  where  the  cooler  spinning  air  attempts  to  penetrate 
the  hot  pilot  exhaust  gases.  This  forces  intimate  contact  between  the  fuel  vapor  and  the  available 
oxygen  at  temperatures  sufficient  for  auto  ignition.  The  rate  of  combustion  is  primarily  driven  by 
the  aerodynamic  forced  mixing  rather  than  locally  turbulent  flame  speed.  Although  the  reaction 
proceeds  in  a fuel  rich  gas,  the  significant  kinetics  are  controlled  by  the  interfacial  composition 
between  the  fuel  vapor  closed  and  the  air  core  of  the  vortices. 

The  subscale  test  data  are  typically  represented  by  Figure  90.  The  system  is  sensitive  to  the 
secondary  air  swirler  angle  since  the  angle  establishes  the  strength  of  the  penetrating  vortex.  The 
system  is,  however,  quite  insensitive  to  the  overall  fuel-air  ratio. 

Similar  tests  were  run  with  a different  form  of  vortex  generator  in  the  fan  duct  stream. 
These  data  are  shown  in  Figure  91.  Although  slightly  more  decrease  in  efficiency  is  observed  with 
increased  fuel-air  ratio,  the  decrease  is  insufficient  to  drive  rumble. 

From  the  results  available  at  this  time,  the  VORBIX  augmentor  should  exhibit  rumble-free 
operation  at  a significantly  higher  altitude  than  the  conventional  system.  The  model  cannot  be 
developed  much  beyond  this  point  without  more  detailed  information  on  the  basic  processes. 
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System  Efficiency  - •/.  System  Efficiency  - % 


The  following  comments  apply  to  a correctly  designed  swirl  on  VORBIX  augmentor: 

• The  controlling  mechanism  on  relative  flame  speed  has  been  decoupled  from 
a direct  fuel-air  ratio  dependence 

• The  piloting  region  has  been  made  dynamically  stifTer  than  a bluff  body 
recirculation  zone 

• The  effect  of  poor  fuel  vaporization  has  been  reduced  or  eliminated. 

All  of  these  factors  generate  reduced  efficiency  response  to  dynamic  variations  in  fuel-air 
ratio,  pressure,  velocity,  or  airflow,  which  in  turn  causes  lower  rumble  gain  factors  and  improved 
stability. 

5.  CONCLUSIONS  AND  RECOMMENDATIONS 

As  a result  of  tbe  exploratoiy  efforts  performed  under  this  study,  the  major  drivers  of  rumble 
instability  in  turbofan  augmentors  have  been  identified.  The  re.sponse  for  a turbofan  augmentor 
of  mixed  flow  configuration  with  conventional  bluff  body  flameholders  bas  been  quantified  and 
the  major  limitations  identified.  A computer  program  has  been  written  which  predicts  the 
combustion  efficiency  of  a conventional  turbofan  augmentor  and  has  been  evaluated  successfully 
against  development  experience  on  a current  engine.  Tbe  possible  improvements  with  advanced 
augmentation  concepts  have  been  explored. 

As  a result  of  the  analytical  and  experimental  efforts  of  this  study  with  regard  to  flame 
stabilization  and  rumble,  the  following  major  conclusions  have  been  reached. 

• Rumble  is  driven  by  oscillations  in  overall  augmentor  efficiency  caused  by  tbe 
response  of  tbe  combustion  process  to  variations  in  fuel-air  ratio,  pressure, 
velocity,  or  temperature.  These  are  important  in  this  order. 

• Consideration  of  the  physical  processes  which  control  flame  stability  in  a two- 
pbase  fuel  spray  bas  resulted  in  the  identification  of  surface  vaporization 
from  the  flameholder  as  a major  parameter  in  the  stability  and  rumble  limits. 

• Tbe  flame  stability  requirements  ultimately  control  the  response  of  the 
overall  augmentation  process  through  the  nature  of  the  turbulent  flame 
initiation  process  in  the  flameholder  wake  shear  layers. 

• Rumble  irstabilities  are  initially  driven  by  a severe  decline  in  tbe  fan  duct 
flameholding  capability  as  the  wake  fuel-air  ratio  approaches  the  rich  limit. 

Due-  to  the  nature  of  the  two-phase  fuel  stability  equations,  this  decline 
occurs  at  total  fuel-air  ratios  well  below  the  rich  limit  overall. 

• Local  modifications  to  the  stability  level  of  the  fan  duct  wakes  through 
flameholder  alterations  or  external  heat  addition  can  greatly  improve  the 
rumble  onset  limits. 
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• Rumble  response  may  be  reduced  by  isolating  the  flame  speed  from  local 
reaction  rate  or  piloting  strength  and  by  improving  the  “hardness”  of  the 
piloting  source.  Two  novel  augmentor  concepts  jifer  the  promise  of  increased 
altitude  capability  without  rumble.  The  VORBIX  offers  the  greater  capabili- 
ty. 

• The  conventional  flameholder  concept  in  a mixed  flow  turbofan  augmentor 
often  limited  improvement  unless  more  control  over  t’!“  wake  composition  is 
obtained  or  methods  are  employed  to  increase  the  wake  reaction  rale. 

The  augmentor  model  as  currently  developed  presents  a reasonable  representation  of  the 
combustion  process  in  the  turbofan  augmentor.  For  an  increased  rumble  limit  conventional 
design,  the  following  guidelines  are  recommended; 

• The  fuel  system  and  flameholder  design  should  be  speciflcally  tailored  to  the 
augmentor  inlet  profile  to  eliminate  locally  rich  fuel-air  regions  and  operate 
at  a uniform  level  of  wake  aerodynamic  loading.  These  profiles  should  he 
realistic  engine  data  as  experience  has  shown  that  these  are  rarely  periodic  or 
symmetric. 

• A high  level  of  wake  reaction  speed  should  be  maintained  by  direct  engine 
stream  mixing  or  heat  addition  to  the  wakes. 

• The  fuel  system  should  be  as  closely  coupled  to  the  wake  as  feasible  and 
consistent  with  acceptable  lean  limits  to  reduce  excessive  wake  enrichment. 

• The  fuel  atomization  should  be  as  great  as  possible  (small  droplets)  for  the 
same  reasons. 

• Methods  of  decoupling  the  flame  speed  from  the  local  fuel-air  ratio  offer  the 
most  significant  areas  of  future  augmentor  design. 

• The  predictions  of  this  model  should  be  verified  on  full-scale  engine 
evaluation. 
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APPENDIX  A 
TEST  DATA 

This  appendix  contains  the  tabulated  and  plotted  experimental  rig  data.  The  following 


symbols  are  used: 

T air  temperature  upstream  of  the  flameholders  ®F 

P rig  pressure  upstream  of  the  flameholders  psia 

M Mach  number  upstream  of  the  flameholder 

ri  combustion  efficiency  ‘x 

TFHl,  TFH2  center  flameholder  metal  temperatures  ®F 


Amp  peak-to-peak  pre.ssure  amplitude  at  flameholders  '"r 

Freq  frequency  read  from  o-graph  Hz 

Included  with  this  tabulation  are  the  transfer  functions  for  the  various  pressure  probes.  The 
following  symbols  are  used  to  define  the  transfer  functions: 

P31.  Po.  P<j.  P«  - ratio  of  pressure  amplitude  at  the  subscript  location  to  the 
amplitude  at  the  flameholder. 

021.  031.  042.  043.  044  * phsse  difference  between  the  pressure  at  the  subscript 
location  and  the  flameholder. 

fai.  fsi.  ^42.  1^43.  1^44  - transfer  function  frequency. 
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TABLE  A-1 

experimental  study  rig  data 


table  A-1 

experimental  study  rig  data  (CONT.) 


EXPERIMENTAL  STUDY  RIG  DATA  (CONT.) 
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EXPERIMENTAL  STUDY  RIG  DATA  (CONT.) 
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EXPERIMENTAL  STUDY  RIG  DATA  (CONT 
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EXPERIMENTAL  STUDY  RIG  DATA  (Cc  il.) 
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EXPERIMENTAL  STUDY  RIG  DATA  (CONT. 
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EXPERIMENTAL  STUDY  RIG  DATA  (CONT.) 
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Figure  A-I.  Bumble  Data  for  Test  Point  Number  t and  2 — Test  Number  7.01 


Configuration 
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•2.  Rumble  Data  for  Test  Point  Numbers  — Test  Number  11. 01 


Average 
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Figure  A-5.  Rumble  Data  fo'  Test  Point  Nuinoer  / — Test  Number  8.01 


Rumble  Data  for  Test  Point  Number  8 — Test  Number  11.01 


Configuration 
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Figure  A-9.  Rumble  Data  for  Test  Point  Number  II  — Test  Number  14.01 


APPENDIX  B 
EQUATIONS 


The  equations  which  are  referred  to  in  the  text  are  presented  here  for  reference. 


Wf,  = diT^f  B-1 

Wfv  = B-2 

Wf,  = jSiWf  B-3 

w,  = p,VA  B-4 

Wf^  = 6 w,  B-5 


Wf,  = = ffiO  w, 

Wf,  = w,.^  - Wf,  = (1-fit)  e w. 
w.  r = rw,  = Fp,  VA 
w,,^=  (l-fii)  « w.r 

Wf,r=  /9i«w.r 


(f/a)* 


Kr 

w.r 


fitO  w.r 
w.r 


= fit« 


(f/a),  = (1-fit)  e 


W,^  = diWf.r  = fitd-fit)  0 w.r 
W = fitVf,,^  = /3,  di(l-/9|)  « w.r 
w = Wf^pK,  = fiiOw,  r K, 

Wf„  = ^i^.rK,  + (i-fit)fi,fi,9w,r 
w.^  = w.rK, 


<f>  = fit0  + (l-fit)0 


Ml 

K. 


-^=  fit  + (1-/J.) 


fiA 

K. 


B-6 

B-7 

B-8 

B-9 

B-10 

B-11 

B-12 

B-13 

B-14 

B-15 

B-16 

B-17 

B-18 

B-19 

B-20 
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EQUATIONS  (ConUnued) 


dl  = dlT  + 

B-21 

Nu„  = 2 + 0.6  Re''>  Pr‘« 

B-22 

Num  = 2 + 0.6  Re''*  SC'* 

B-23 

w = KA.p.ln  (p^P  p ) 

B-24 

„ Nu  Dv  MW 

R4T. 

B-25 

L k Nuh 

’’’■  i, 

B-26 

q = h,A,(T.-T, ) /3 

B-27 

^ = e*-l 

B-28 

2 = Cpv  w/irkd(  Nuh 

B-29 

Aq  = q - w»X 

B-30 

dX  Aq  Aq 

dt  ni/  Cpi  4/3  *•  (d(/2)*  p,  Cp»  ‘ 

(31) 

^ = 4-^^  <v.-v,.- 

B-32 

„ P.d,(V.-V,) 

Re  

B-33 

Wf,  = w,j  rdi 

B-34 

T*r 

B-35 

Wfc  = v, 

B-36 

w,  = pV,  A 

- V N 

~ Pl^l  IT 

A - N 

B-37 
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EQUATIONS  (Continued) 


10 

= - — A H m,  , Afj„ 
m,  , A , , 


w„ 

Wv 

Pv 

T 

q = 

W, 


dm 
dt 

= fcn(Num,  Pv.  A.) 

= fcn(T,  ) 

= T, , + ‘ i at. 

^ w,X  + w,,.  C,.  AT. 

^ fen  (Nuh,  T„  T*) 

■ s^) 


C.  = 


pH  Pv 

Nun,  Dv  MW 


Nuw  = 


R AxT. 
h N 

-r— = 0.99  Re"  * P,"  ** 

Km 


q = 


w,,  C,  ,T  ,-T  A.p.l„  ( ) 


W.,  = P.V.NK, 


V„  = 


= ilsXi 

T 

^ P.  V„  /r 
Cv(L/D)(B/D)N« 
rV. 

" ^-W 

w = 

N r- 

p,V.Cv(L/D)(B/D)N 

'-r 

K,  = Cv{L/D)(B/D)(r')-> 


B..38 

B-.39 

B-40 

B-41 

B-42 

B-4.3 

B-44 

B-45 

B-46 

B.47 

B-48 

B-49 

B-50 

B-5I 

B-.52 

B-53 

B-54 

B-55 

B-lie 
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EQUATIONS  (Continued) 


P.NT,'  • ' 

p.V.C,(L/D)(B/D)N 


A = w.^ 


V„  = C.(UD)(B/D)N» 


- JEl_ 

RT. 


A 

VoP* 

A 

V„p* 

A 

VoP* 

dmo 


RT.  r'N  p.* 


1 Va 

tR  P.N 


= constant  X 


V. 


P.NT.‘ 


k e'®'''' 

dt  “ "rT^  *"*  *'"  * 

A _ k(ni+l) 


VoP" 


R."y«  X.‘ Xf’" 


0CtH(  + 50j  + 5niN|  — • Syi  CO*  + 4y«H,0  + 

(0-y«)C,H,  + 5(l-y«)  0,  + 5mN, 

T = T,  + «AT 

A _ k(m+l)  [5(l-y«)l*  [0-y«l"~" 

VoP"  R,"y«  [5{m+l)  + 0 + ytj"  (T,+«AT1"-** 

A ^ 1.29  X 10‘«  (m+l)  [5(l-y€)l*U-yt)*e-®^T,>.AT) 
V„p**  (0.08206) ’*ytf5(m+l)  + ^ + y€j»*IT,+«ATl 


VoP" 

A 


430  k, 


e-«'T,  {20(l-«)r  (l-0«) 


T,"  ^[4.76  + 0(1.36- *)]>• 

430  . (0.080  0)»«/ 1-t  1" 

« K*  j ,,t  ^ I 4.76-*  + 0.080  (1  + 16«)' 


VoP 

(-^) 

V VoP*  / mu.  pnd 

' VoP*  ' BUI,  dBU 


= 3.55 


B-57 

B-58 

B-59 

B-60 

B-61 

B-62 

B-63 

B-64 

B-65 

B-66 

B-67 

B-68 

B-69 

B-70 

B-71 

B-72 
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EQUATIONS  (Continued) 


U = 


V. 


(l-D 


St  = + (2uS> )" 

The  Value  of  u'  is: 


B-73 

B-74 

B-75 


u'  = «,U 

B-76 

Sf  = St  • I,. 

B-77 

II 

o 

n-78 

AYo  = If,,  • w 

B-79 

i^y 

B-80 

B-81 

-J-=  (0.32)  + (0.68) = 1.325 

B-82 

-1-=  d.  + (1-^.)-^ 

B-8.3 

53-  angle:  (0.32)  + (0.68) = 1.325 

B-84 

90»  angle:  |-=  (0.32)  + (0.68)  = 0.912 

B-85 

= d,  + (1-/3.)  ■ 

B-86 
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EQUATIONS  (Continued) 


V n«  “ 
oH 

B.87 

0 = 00,  + 0(1-0,) 

B-88 

d,  = fen  (T..  0„  e) 

B-89 

Tw  “ T,  + ij, 

B-90 

= fen  (T.,  0) 

B-91 

tiw  = fen  (T„  A/V„p‘) 

B-92 
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LIST  OF  SYMBOLS 


■ \ 


English 

a Reac*ion  order  for  oxygen  (eqn.  64  > 

A Duct  flow  area 

A Stirred  reactor  maM  loading,  gm-moie/sec 

A,  Surface  area  for  heat  or  mass  flux 

B/D  Wake  width/flameholder  width 

C Kinetic  activation  constant  = E/R 

Ctf  Drag  coefficient 

C,  Wake  shape  factor  (eqn.  52) 

Cp  Specific  heat 

C.  Surface  mass  efflux  parameter  (eqn.  46) 

Dv  Binary  diffusion  coeflicient 

4 Liquid  fuel  droplet  diameter 

E Activation  energy 

f/a  Fuel-air  ratio 

hr  Film  coefficient  for  heat  flux 

k Film  thermal  conductivity 

K Mass  efflux  parameter  (eqn.  25) 

k|.  k*  Frequency  factor  for  kinetic  collision  theory 
Wake  recirculatioo  coeflicient  (eqn.  17) 

L/D  Wake  lengthfllameholder  width 

m Ratio  of  diluent  to  oxygen  mole  fraction 

m Maas  quantity 

in«  Mass  of  oxygen  in  wake 

MW  Molecular  weight 
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LIST  OF  SYMBOLS  (Continued) 

n Overall  reaction  order 

N Flameholder  width 

Nuh  Nusselt  number  for  heat  flux 

Num  Nusselt  number  for  mass  flux 

Pr  Prandtl  number 

P Pressure 

q Heat  flux 

R,  R,  Gas  constant 

Re  Reynold’s  number 

Sc  Schmidt  number 

S Laminar  flame  speed 

St  Turbulent  flame  speed 

t Time 

T Temperature 

Tr  Flame  temperature  (eqn.  70) 

u RMS  turbulence  velocity 

U Accelerated  air  velocity  (eqn.  73)  " 

V Velocity 

Vo  Reactor  volume 

w Flowrate 

W Duct  width 

X Axial  distance 

y Reactor  flow  factor;  = ^ if  ^ < 1 and  = 1 if  ^ > 1 

z Blocking  coeflicient  (eqn.  29) 


LIST  OF  SYMBOLS  (ContinuMf) 


Greek 

a Flamehnldrr  aprx  anjslr 

di  pTF-flamaholder  vapnti^atinn  ctiefTirirnt  (rqn.  3) 

dt  F)irl  odlaction  mrffirH>nl  <aqn.  13l 

di  Surface  vap«mrati«>n  ctiefTirirni  (acjn.  14> 

d,T  Thnittt'.ni;  pmcewt  X'aponzalinn  rtwfTicir'it  <eqn.  1 1 

d,,  Dmplel  vaporizatMtn  rttefTirient  leqn.  2> 

d Blorkinj!  oieflicient  (FC)n.  2St 

r BInckacF  ratio 

*•  Capture  atreamline  ieqn.  3(;> 

# Wake  nzvfrm  cnnauir.pt  ion  efficiencx- 

« Turbulence  intenaiiy  at  flameholder  lip 

f Wake  reaction  efficiency 

# Overall  fuel -air  ratio 

A Latent  heat  uf  vaporization 

0 NlacrMity 

X.  Oiypen  tnolal  concentration 

Xt  Fuel  molal  concentration 

# Wake  fuel-air  ratio 

r Residence  time 

r Non-dimensional  residence  time 

a Deitaity 


iS2 


LIST  OF  SYi'.ioCLS  (Continued) 


Su  hacripts 
a Air 

f Fuel 

' Liquid  phase  fuel 

^apor  phase  fuel:  vaporized  portion 
t Total 

c Collected  on  flamcholder 

r Pertaining  to  flanieholder  area 

w Into  the  wake  of  the  flameholder 

s Static,  surface,  sensible 

0 Initial  or  injection  conditions 

1 Initial  or  stepping  notation 

r Recirculated 


Operational  Notes 

A Finite  difference 

d Differential  operation 

a Partial  differential  operation 

- Average  value 
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